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Volume  II 


VALIDATION  OF  THERMOCHEMkCAL  ANALYSIS  TECHNIQUES 
FOR  HIGH  CHAMBER  PRESSURE  SOLID  ROCKET  NOZZLES 
WITH  HIGH  VELOCHY  BLAST  TUBES 


by 


D.  L.  Baker 
A.  J.  Murphy 
M.  R.  Wool 


The  work  discussed  in  volume  II  -  was  performed  by 
Aerotherm  Division  of  Acurex  Corporation  under 
subcontract  to  United  Technology  Center  on  contract  No. 
F0461 1-71-C-0051.  This  report  presents  the  final  work  of  the 
project. 
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ABSTRACT 


The  Jet  of  Aerotherm  computer  programs  (ACE,  CMA,  ARGEIBL)  which  were 
validated  ior  charring  ablating  and  carbon/carbon  materials  in  conventional 
-~nd  low  velocity  blast  tube  nozzles  under  Contract  F04611-69-0-0065’  were 
validated  for  these  same  materials  for 

•  high  velocity  blast  tube  nozzles 

•  operation  at  high  chamber  pressures  (up  to  3000  psia) 

•  aluminized  CTPB  propellants  which  vary  in  aluminum  content  from  5  to 
18  percent. 

The  validation  was  determined  based  on  comparisons  between  measured  and 
predicted  surface  recession  and  char  depth  thicknesses  and  between  predicted 
and  determined  convective  heat  transfer  coefficients,  ^he  GASKET  computer  program 
which  considers  kinetically  controlled  surface  thermochemical  reactions  was  used 
for  predicting  the  performance  of  the  bulk  graphites  and  pyrolytic  graphite 
materials.  This  codf  was  validated  for  the  bulk  graphites  but  not  for  the 
pyrolytic  graphite  washers  in  this  particular  environment.  Under  Contract 
F04611-69-C-0081,  this  code  had  been  validated  for  both  pyrolytic  graphite  washers 
and  coatings  in  a  wide  variety  of  environments. 

The  surface  recession  and  char  depth  measuremer,  .s  used  for  comparisons  with 
the  analytical  predictions  were  obtained  from  materials  located  in  the  aft 
closure,  blast  tube,  aft  entrance  cap,  nozzle  throat,  and  exit  cone  of  nozzles 
which  were  exposed  to  a  chamber  pressure  of  approximately  3000  psia  for  10  seconds. 
This  analysis  volume  summarizes  the  initial  13  test  firings  results.  Three 
additional  tests  were  conducted  at  chamber  pressures  up  to  3500  psia  using  wire 
wound  tungsten  inserts.  The  results  for  these  last  three  tests  are  presented  in 
Volume  I  of  this  report.  The  rocket  motor  used  for  these  tests  nad  either  a 
cylindrical  or  keyhole  propellant  grain  configuration  with  either  a  5  percent 
Al-CTPB  or  18  percent  Al-CTPB  propellant.  The  convective  heat  transfer  coefficient 
used  for  comparison  with  the  analysis  techniques  was  obtained  from  two  total  heat 
flux  sensors  installed  in  the  blast  cube  of  one  of  the  rocket  motor  firings. 

In  general,  the  agreement  between  the  analytical  predictions  and  the  surface 
recession  and  char  depth  data  was  either  within  t  25  percent  or  within  the  ex¬ 
perimental  accuracy  of  the  measurements  with  the  exception  of  the  throat.  The 
best  agreement  was  obtained  for  the  analysis  locations  in  the  blast  tube  and  exit 
cone.  The  poorest  agreement  occured  for  the  pyrolytic  graphite  washers  in  the 
nozzle  throat  region. 
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SECTION  1 


INTRODUCTION 

The  wide  variety  of  mission  requirements  for  high  chamber  pressure  solid 
rocket  motors  result  in  a  correspondingly  wide  range  in  the  operational  character¬ 
istics  of  the  motor  and  nozzle.  This  wide  range  in  the  motor  and  nozzle  charac¬ 
teristics  dictates  that  a  relatively  broad  material  performance  technology  base 
be  available  in  order  for  the  engineer  to  efficiently  design  a  nozzle  which  will 
satisfy  the  stringent  thermal  and  structural  requirements  for  a  flightweight 
nozzle  design.  The  motor  performance  requirements  of  particular  interest  to 
the  nozzle  design  engineer  are  the  motor  ballistics  and  the  type  of  propellant. 
Depending  on  the  particular  mission  requirement,  the  high  chamber  pressure  motor 
ballistics  can  range  from  a  short  to  moderate  single  pulse  motor  operation  to 
a  multiple  pulse  or  boost  sustain  type  of  motor  operation.  An  aluminized  pro¬ 
pellant  is  generally  used  in  these  motors  but  the  aluminium  content  in  the  pro¬ 
pellant  can  vary  from  5  to  18  percent.  Typical  nozzle  operating  characteristics 
can  include  restrictions  on  the  nozzle  configuration  as  well  as  on  the  nozzle 
performance.  The  nozzle  configuration  restrictions  may  require  a  high  velocity 
blast  tube  and  the  nozzle  performance  restrictions  generally  include  surface 
recession  limits  (particularly  at  the  nozzle  throat)  and  limitations  on  the 
temperature  rise  of  the  nozzle  structural  shell. 

Because  of  the  broad  material  performance  technology  which  is  required 
to  satisfy  the  wide  range  of  material  and  nozzle  performance  requirements,  several 
Air  Force  programs  have  been  conducted  with  the  objectives  of  1)  obtaining 
experimental  material  performance  data  at  high  chamber  pressures  and  2)  using  these 
data  to  validate  existing  analytical  techniques.  These  programs  included  the  mat¬ 
erial  performance  evaluation  of  carbon  phenolic  ablators  and  pyrolytic  graphite 
(PG)  coatings  when  subjected  to  multiple  pulse  duty  cycles  under  Contract  F04611- 
70-C-0019  (Ref.  1),  the  definition  and  evaluation  of  the  kinetic  response  of  pyro¬ 
lytic  graphite  under  Contract  F04611-69-C-0081  (Ref.  2)  and  the  material 
performance  evaluation  of  several  charring  ablators,  carbon-carbons,  bulk  graphites, 

PG  washers  and  wire  wound  tungsten  under  Contract  F04611-69-C-0065  (Refs.  3  and  4). 
Under  all  of  wiese  contracts,  a  set  of  analytical  techniques  were  validated  based 
on  comparisons  between  measured  and  predicted  material  performance.  The  objectives 
of  this  contract  (F04611-71-C-0051)  are  to  provide  material  performance  data  and  to 
update  the  validation  of  the  analytical  techniques  for  a  wide  range  of  aluminum 
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loadings  in  the  solid  propellant  {5  to  18  percent)  and  for  high  as  well  as  low 
velocity  blast  tubes  (blast  tube  area  ratio  from  1.3  to  2.0  with  respective  Hach 
numbers  of  .53  and  .31. 

In  support  of  Contract  F04611-71-C-0051  and  under  subcontract  to  UTC, 
Aerotherm  performed  thermochemical  and  thermal  analyses  and  provided  thermal  in¬ 
strumentation.  The  specific  items  covered  in  this  subcontract  were 

Phase  I  Thermochemical  screening  calculations  of  ablation  materials 
to  determine  their  performance  in  a  high  chamber  pressure 
solid  rocket  motor  environment. 

In-depth  thermal  response  predictions  using  the  results  of 
the  thermochemical  screening  calculations. 

Phase  II  Definition,  design,  and  fabrication  of  thermal  instrumentation 
for  installation  into  a  high  chamber  pressure  motor. 

Transient  two-dimensional  heat  conduction  analyses  of  a 
pyrolytic  graphite  washer  nozzle  throat  package. 

Phase  III  Eighteen  material  performance  predictions  which  include 

surface  recession,  char  depth  regression  (if  applicable) , 
and  in-depth  transient  temperature  histories. 

Evaluation  of  the  method  for  predictina  the  convective  heat 
transfer  coefficient. 

Phase  IV  Post-test  experimental  evaluation  of  the  material  density 
profile. 

The  analytical  techniques  used  in  Phases  I  thorugh  III  are  summarized 
in  Section  2.  The  material  property  data,  propellant  property  data,  nozzle 
configuration,  and  the  convective  and  radiative  heat  flux  boundary  conditions 
are  presented  in  Section  3.  The  results  obtained  for  the  thermochemical  screen¬ 
ing  and  in-depth  thermal  response  predictions  and  for  the  pyrolytic  graphite 
washer  throat  package  are  presented  in  Sections  4  ana  5,  respectively.  The 
analytical  material  performance  predictions  for  the  motor  firings  are  presented 
along  with  comparisons  with  the  experimental  data  in  Section  6.  The  measured 
results  from  the  thermal  instrument ation  are  reported  in  Section  7  and  finally, 
conclusions  and  recommendations  ara  presented  in  Section  8. 
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SECTION  2 


SUMMARY  OF  ANALYSIS  TECHNIQUES 

As  mentioned  previously  the  objectives  of  three  previous  Air  Force 
Contracts  (F04611-69-C-0065 ,  F04611-70-C-0019  and  F04611-69-C-0081)  included 
the  validation  of  analytical  techniques  for  predicting  the  thermochemical  res¬ 
ponse  of  materials  which  are  typically  used  in  high  chamber  pressure  solid 
rocket  motors.  Under  Contract  F04611-69-C-0065,  the  comparisons  between  predicted 
ani  measured  surface  recession  showed  that  the  analytical  techniques  predicted 
the  surface  recession  of  charring  ablative  materials,  carbon-carbon  materials, 
and  wire  wound  tungsten  within  t  25  percent.  Under  Contract  F04611-69-C-0081, 
similar  comparisons  showed  that  when  the  analytical  techniques  were  modified 
to  include  kinetically  controlled  surface  chemical  reactions  and  if  the  proper 
kinetically  controlled  constants  were  available,  then  che  surface  recession 
histories  for  PG  coatings  and  PG  washers  were  generally  predicted  within  ±  50 
percent  for  a  wide  variety  of  thermochemical  environments.  These  two  contracts 
considered  only  single  pulse  motors.  Under  Contract  F04611-70-C-0019 ,  compari¬ 
sons  between  the  predictions  and  measurements  showed  that  equally  good  results 
were  obtained  from  multiple  pulse  duty  cycle  motors  for  carbon  phenolic  and  PG 
coatings.  Detailed  results  of  these  studies  are  presented  in  References  1 
through  4.  For  the  thermochemical  analyses  performed  under  the  current  contract, 
these  same  analytical  techniques  were  used,  and  they  are  summarized  briefly  in 
the  following  paragraphs. 


The  analytical  prediction  techniques  are  in  the  form  of  computer  programs 
and  are  identified  by 

•  Aerotherm  Chemical  Equilibrium  (ACE)  Program 

•  Charring  Material  Ablation  (CMA)  Program 

•  Axisymmetric  Transient  Heat  and  Material  Ablation  (ASTHMA)  Program 

•  Aerotherm  Real  Gas  Energy  Integral  Boundary  Layer  (ARGEIBL)  Program 

•  Aerotherm  Graphite  Surface  Kinetics  (GASKET)  Program 

•  Boundary  Layer  Inversion  Matrix  Procedure  (BLIMP)  Program 

•  Steady  State  Ablation  (SSA)  Program 
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The  SSA  Program  which  predicts  the  steady  state  recession  and  surface  tempera¬ 
ture  of  a  material  was  used  along  with  the  ACE  and  ARGEIBL  programs  for  making 
the  thermochemical  screening  predictions  presented  in  Section  4.  The  analysis 
procedure  for  making  these  predictions  is  summarizad  in  Figure  2-1.  The  first 
five  programs  were  the  basic  ones  used  in  making  the  analytical  predictions  which 
are  presented  in  Sections  5  and  6.  The  analysis  procedure  for  making  these 
predictions  is  summarized  in  Figure  2-2.  As  described  in  Reference  3,  the 
BLIMP  program  is  used  primarily  to  evaluate  the  heat  and  mass  coefficients 
predicted  by  the  more  simplified  analytical  techniques.  In  this  reference,  the 
BLIMP  proaram  was  used  to  confirm  the  heat  and  mass  transfer  coefficients  pre¬ 
dicted  by  the  much  simpler  ARGEIBL  program  and  Chilton-Colburn  correlation  for 
a  typical  high  chamber  pressure  nozzle  configuration.  The  major  features  of 
the  ACE,  CMA,  ASTHMA,  ARGEIBL,  GASKET,  and  BLIMP  computer  programs  are  summarized 
in  Table  2-1.  For  a  more  detailed  description  of  the  programs,  the  reader  is 
referred  to  References  4  through  8. 

Versions  of  the  ACE,  CMA,  and  ASTHMA  computer  codes  which  were  used  in 
making  the  analytical  predictions  are  available  to  the  rocket  community.  The 
version  of  ACE  which  is  available  is  the  EST3  code.  This  code,  which  is  described 
in  greater  detail  ir.  Reference  9,  treats  the  diffusion-controlled  thermochemicax 
response  of  materials,  and,  therefore,  it  is  applicable  to  all  materials  for 
which  the  surface  chemical  reactions  are  not  kinetically  controlled  or  for 
which  a  surface  melt  layer  does  not  develop.  These  materials  include 
carbon  and  graphite  phenolic  charring  ablators  and  carbon/carbon  composites. 
However,  EST3  is  not  applicable  for  graphites  and  silica  phenolic  charring 
ablators.  The  CMA  computer  code  described  in  Table  2-1  is  available  as  CMA3 
(Reference  10) ,  The  ASTHMA  and  GASKET  codes  described  in  Table  2-1  were  made 
available  through  Contract  F04611-69-0081 ,  and  the  background  to  these  codes 
are  presented  in  References  11  and  12. 


TABLE  2-1 


SUMMARY  OF  24AJOR  FEATURES  OF  COMPUTER  CODES 


Computer 

Program  Major  Features  Typical  Use 

CMA3  1  Calculation  of  transient  thermal  1.  Thermal  analysis  of  ablative 
response  of  a  composite  slab  or  materials  in  liquid  and  solid 

cylinder  containing  up  to  five  in-  rocket  motors, 

depth  charring  (or  noncharring) 
materials . 

2.  Calculation  of  surface  recession  2.  Thermal  analysis  of  reentry  nose 

rate  resulting  from  diffusion  or  tips. 

kinetics.lly  controlled  surface 
chemical  reactions  and/or  mechan¬ 
ical  melt  removal. 

3.  Considers  variable  thermal  proper¬ 
ties  with  a  flexible  physical  model  of 
thermal  conductivity  which  permits 
the  treatment  of  partially  charred 
values  less  than  the  virgin  material 
value. 

4.  Considers  in-depth  and  surface 
thermochemical  effects  of  the 
pyrolysis  gases  being  generated  by 
both  the  surface  and  backup 
materials. 


GASKET  1.  Calculation  of  thermodynamic  state  1.  Nozzle  equilibrium  isentropic 
for  isentropic  or  isoenthalpic  expansion, 

process  for  a  closed  system 
including  condensed  phase. 

2.  Calculation  of  transport  proper-  2.  Evaluation  of  transport  coeffici- 

ties  (viscosity,  thermal  conduc-  ents  for  use  in  boundary  layer 

tivity,  and  mass  diffusion  analysis. 

coefficients) . 

3.  Calculation  of  surface  state  of  3,  Generation  of  Mollier  charts  for  any 

graphite  materials  exposed  to  an  gas  system. 

environment  considering  kinetically 
controlled  reactions  between  the 
material  and  the  H,0,  CO,,  and  H2 
of  the  combustion  products.  “ 

4.  Considers  unequal  diffusion  4.  Generation  of  thermochemical  data 

coefficients.  for  use  with  CMA  and  ASTHMA  computer 

programs 
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TABLE  2-1  (Concluded) 


ASTHMA 3  1.  Calculation  of  transient  thermal  1.  Thermal  analysis  of  heat  sink 

response  of  a  two-dimensional  materials  in  liquid  and  solid  rocket 

body.  motors  and  reentry  nose  tips. 

2.  Considers  surface  recession 

3.  General  thermochemical  surface 
bovudary  conditions 

4.  Co;  aiders  anisotropic  temperature 
dependent  thermal  properties. 


ACE  1.  Cilculation  of  thermodynamic  1.  Nozzle  equilibrium  isentropic 

state  for  isentropic  or  isoenthal-  expansion, 

pic  process  for  a  closed  system 
including  condensed  phases. 

2.  Calculation  of  transport  properties  2.  Evaluation  of  transport  coefficients 

(viscosity,  thermal  conductivity,  for  use  in  boundary  layer  analysis, 

and  mass  diffusion  coefficients). 

3.  Calculation  of  surface  state  of  any  3.  Generation  of  Mollier  charts  for  any  gas 

material  exposed  to  an  environment  system. 

considering  all  possible  chemical 
reactions  between  the  materials  and 
the  environment.  Surface  reactions 
may  be  kinetically  controlled. 

Mechanical  removal  of  surface  melt 
layer  may  also  be  considered. 

4.  Considers  unequal  diffusion  coeffi-  4.  Generation  of  thermochemical  data  for 

cients.  use  with  CMA  and  ASTHMA  computer 

programs. 

ARGEIBL  1.  Calculation  of  laminar  or  turbulent  1.  Transfer  coefficient  analysis  of 
convective  heat  transfer  coeffici-  solid  and  liquid  rocket  motors, 

ents  for  non-ablating  wall. 

2.  Considers  any  real  gas  system  and 
axial  variations  in  wall  tempera¬ 
ture  and  free  stream  properties. 

3.  Solves  the  energy  integral  equs  .on. 

BLIMP  1.  Calculation  of  laminar  or  turbulent  1.  Transfer  coefficient  analysis  of 

convective  heat  transfer  coeffi-  solid  and  liquid  rocket  motors  and 

cients  for  ablating  or  nonablating  reentry  nose  cones. 

wall. 

2.  Considers  any  real  gas  system 

and  axial  variations  in  wall  tempera¬ 
ture,  free  stream  properties,  and 
surface  recession. 

3.  Considers  chemical  and  thermal 
diffusion  in  the  boundary 
layer . 

4.  Solves  nonsimilar  laminar  and  turbulent 
boundary  layer  equations. 
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FIGURE  2-1.  FLOW  CHART  FOR  SCREENING  ANALYSIS 


.  ACE  OR  GASKET  CODE 


Fiaure  2-2.  Flow  Chart  for  Detailed  Transient  Analysis 


SECTION  3 


PRELIMINARY  ANALYSIS 

In  the  previous  section,  the  analysis  techniques  used  to  study  the  response 
of  ablative  rocket  nozzles  to  high  pressure  combustion  products  were  presented. 

The  first  step  in  the  analysis  procedure  consists  of  defining  the  free  stream 
properties  by  making  computations  using  the  ACE  computer  code.  The  next  step  con¬ 
sists  of  defining  the  convective  heat  transfer  coefficients  by  making  computations 
using  the  ARGEIBL  code.  The  results  of  these  calculations  which  depend  only  on 
the  propellant  properties  and  nozzle  geometry  are  presented  in  this  section. 

They  were  compiled  for  the  two  nozzle  geometries  (blast  tube  area  ratios  of  1.3 
and  2.0)  and  the  two  propellants  (5  and  18  percent  aluminum  contents)  used  in  the 
test  program  Dy  UTC.  In  later  sections,  more  specific  calculations  are  pre¬ 
sented  which  define  the  material  response,  and  these  calculations  use  the  re¬ 
sults  of  this  section  for  boundary  conditions  and  material  properties. 

Propellant  data  and  basic  nozzle  geometries  are  presented  in  Section  3.1 
and  3.2,  respectively.  Section  3.3  describes  the  flow  field  analysis  used  for 
the  cylindrical  and  keyhole  port  propellant  cartridges.  The  free  stream 
properties  and  heat  transfer  boundary  conditions  are  presented  in  Sections  3.4 
and  3.5,  respectively.  The  material  property  data  are  defined  in  Section  3.6. 

3.1  FROPELLANT  DATA 

Two  variations  of  a  CTPB  propellant  were  used  for  the  following  analyses. 
One  had  a  5  percent  aluminum  loading  fraction,  and  the  other  had  an  18  percent 
aluminum  loading  fraction.  The  chemical  elements  in  the  5%  Al-CTPB  are 

Element  Atomic  »?ts/Unit  Mass 

Hydrogen  0.0412 

Carbon  0.0083 

nitrogen  0.0071 

Oxygen  0.0287 

Aluminum  0.0019 

Chlorine  0.0071 

The  flame  temperature  for  this  propellant  is  5816PR  at  1000  psia.  The  chemical 
elements  constituents  in  the  18%  Al-CTPB  propellant  are 

Element  Atomic  Wts/Unit  Mass 

Hydrogen  0.0368 

Carbon  0.0082 

Nitrogen  0.0060 

Oxygen  0.0242 

Aluminum  0.0066 

Chlorine  0.0060 

The  flame  temperature  for  this  propellant  is  6469°R  at  1000  psia. 
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3.2 


NOZZLE  GEOMETRY 


The  calculation  of  the  heat  transfer  boundary  conditions  requires  the  de- 
finition  of  the  internal  nozzle  contour  and  a  nominal  expected  surface  tempera¬ 
ture  range  for  the  surface  materials.  Figure  3-1  illustrates  the  nozzle  geome¬ 
tries  which  were  used  for  the  analyses  presented  in  Sections  4,  5,  and  6.  As 
shown  by  this  figure,  two  basic  blast  tube  contraction  ratios  were  considered. 
The  nozzle  geometries  are  identical  for  the  aft  closure,  the  nozzle  throat  dia¬ 
meter,  and  the  exit  cone.  The  primary  difference  in  the  geometries  was  that  one 
had  a  blast  tube  area  ratio  of  1.3  and  the  other  one  had  a  blast  tube  area  ratio 
of  2.0.  This  aifferer.ee  in  the  area  ratio  required  3ome  minor  variations  in  the 
forward  and  aft  entrance  caps  and  in  the  convergence  angle  of  the  PG  washers 
or  wire  wound  tungsten  upstream  of  the  nozzle  throat. 

The  materials  considered  for  each  of  the  nozzle  components  are  described 
in  Section  4  for  the  steady  state  thermochemical  screening  analyses,  in  Section 

5  for  the  two  dimensional  analyses  of  the  nozzle  throat  package,  and  in  Section 

6  for  the  post  test  thermochemicai  predictions.  The  surface  temperature  range 
of  primary  interest  for  these  materials  is  between  2700  and  5490°R. 

It  is  worth  noting  that  the  nozzle  shown  in  Figure  3-1  was  designed  so 
that  components  could  have  more  than  one  material  exposed  to  the  combustion 
products.  The  aft  closure,  for  example,  was  divided  into  two  sections,  one  en¬ 
compassing  the  plane  above  the  centerline  of  the  nozzle  and  the  other  the 
plane  below  the  nozzle  centerline.  Thus,  for  a  single  firing,  two  aft 
closure  materials  can  be  tested.  In  a  similar  manner,  the  aft  entrance 
cap  and  the  exit  cone  had  the  capability  of  testing  two  materials  for  any 
one  firing,  and  the  blast  tube  could  contain  up  to  three  test  materials. 


3.3  FLOW  FIELD  ANALYSIS 

The  purpose  of  the  flow  field  analysis  was  to  define  the  mass  flux  as 
a  function  of  boundary  layer  running  length.  This  allows  the  free  stream  pro¬ 
perties  as  calculated  by  the  ACE  computer  program  to  be  converted  from  a  function 
of  mass  flux  to  a  function  of  nozzle  location.  Two  configurations  of  the  pro¬ 
pellant  loaded  cartridge  were  fired  during  the  course  of  the  orogram.  One 
of  these  configurations  had  a  circular  port  with  burning  both  on  the  aft  end 
and  in  the  port  and  the  other  configuration  had  a  small  circular  port  with  burn¬ 
ing  both  in  the  port  and  in  a  slot  running  the  length  of  the  grain  (Figure  3-2) . 

The  flow  field  for  the  first  configuration  was  analyzed  assuming  one 
dimensional  flow  and  utilizing  the  following  relationship  which  describes  the 
conservation  of  mass  for  internal  flows 
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(3-1) 


-v,  mu  »* 

Pu  =  — tttt  =  r  =  inverse  of  contraction  ratio 
p*Uw  A 
e  e 

where  Peue  and  A  are  the  mass  flux  and  cross  sectional  area,  respectively. 

An  asterisk  represents  nozzle  throat  or  sonic  conditions.  The  relationship 
between  contraction  ratio  and  boundary  layer  running  length  is  shown  in 
Figure  3-3.  By  combi* ing  this  relationship  with  Equation  3-1,  the  variation 
of  the  non-dimensional  flow  parameter  pu  with  the  boundary  layer  running 
length  was  obtained,  aid  this  is  shown  by  the  t  =  “  line  in  Figire  3-4.  The 
slotted  grain  con figur.it. ion,  however,  cannot  be  analyzed  using  one-dimensional 
flow  theory.  The  flow  field  analysis  for  this  configuration  is  described  in 
the  following  paragraphs. 

The  flow  field  for  the  slotted  propellant  configuration  shown  in 
Figure  3-2  changes  rapidly  with  time  due  to  the  increasing  diame<  r  of  the 
circular  port  and  the  widening  of  the. slot.  One  should  note  that  the  propel¬ 
lant  is  inhibited  on  the  aft  face,  and,  thus,  the  surface  of  the  aft  face  of 
the  propellant  does  not  recede  from  the  surface  of  the  aft  closure.  Of 
course,  the  aft  closure  gradually  becomes  exposed  during  the  motor  firing  due 
to  the  opening  of  the  circular  port  and  the  widening  of  the  slot. 

One  of  the  difficulties  in  performing  the  flow  field  analysis  of  the 
slotted  grain  is  in  defining  the  mass  flux  exchange  between  the  mass  of  combus¬ 
tion  products  generated  in  the  slot  with  those  generated  in  the  circular  port. 
At  motor  ignition  the  surface  area  for  the  circular  port  is  smaller  than 
for  the  slot.  Thus,  as  the  flow  of  combustion  products  proceeds  from  the 
forward  to  the  aft  end  of  the  propellant,  the  combustion  products  from  the 
slot  will  feed  into  the  circular  port.  At  later  motor  burn  times,  the 
process  is  reversed;  that  is  the  combustion  products  generated  at  the 
surface  of  the  circular  port  will  feed  into  the  slot.  The  second  difficulty 
in  performing  the  flow  analysis  is  in  defining  the  effective  flow  area  of  the 
slot  as  the  combustion  products  in  the  slot  are  forced  into  the  circular  port 
at  the  aft  end  of  the  motor. 

For  the  analysis  presented  herein,  the  first  difficulty  mentioned  above 
was  resolved  by  making  the  conservative  assumption  that  up  to  the  time  when 
the  combustion  gases  in  the  port  starts  to  feed  the  slot,  all  the  combustion 
products  generated  in  the  slot  flow  down  the  slot  and  that  all  the  combustion 
products  in  the  circular  port  flow  through  the  circular  port.  To  resolve  the 
second  difficulty,  an  engineering  judgement  was  made  in  order  to  define  the 


following  paragraphs  for  the  cylindrical  grain  and  keyhole  port  grains, 
respectively.  For  this  analysis,  a  chamber  pressure  of  2800  psia  was  used. 

The  heat  transfer  coefficient  distribution  in  nozzles  fired  with  a 
cylindrical  grain  propellant  were  calculated  accounting  for  the  effects  of 
1)  both  5  percent  and  1«  percent  aluminum  propellant,  2)  blast  tube  area 
ratios  of  1.3  and  2.0,  3)  and  wall  temperatures  of  2700°  and  5400°R.  Figure 
3-8  illustrates  the  heat  transfer  coefficient  distribution  predicted  by  the  ARGEIBL 
code  for  the  5  percent  aluminum  propellant  for  blast  tube  area  ratios  of  1.3  and 
2.0.  In  each  figure  the  sensitivity  to  wall  temperature  is  shown.  Figure  3-9 
illustrates  the  heat  transfer  coefficient  for  the  18  percent  aluminum  propellant. 
Note  that  the  18  percent  aluminum  CTPB  yields,  in  general,  a  lower  heat  transfer 
coefficient  than  the  5  percent  propellant.  Also,  the  sensitivities  of  each  propel¬ 
lant  to  wall  temperature  are  similar.  For  the  case  of  the  1.3  blast  tube  area 
ratio,  an  analysis  was  performed  to  determine  the  sensitivity  of  the  heat  transfer 
coefficient  on  the  two  dimensional  flow  field  effects  in  the  aft  closure,  and  this 
analysis  is  described  in  the  following  paragraph. 

As  mentioned  previously  (Section  3.3),  a  one  dimensional  flow  field  analysis 
was  used  for  calculations  involving  the  cylindrical  port  propellant  cartridge.  To 
account  for  the  sharp  convergence  angle  of  the  aft  closure,  an  analysis  was 
performed  in  which  the  effective  flow  area  in  the  aft  closure  was  assumed  to 
be  a  spherical  surface  whose  center  is  the  apex  of  the  cone  created  by 
extending  the  aft  closure  surface  to  the  motor  centerline.  This  surface  is 
defined  by  the  schematic  drawing  presented  in  Figure  3-10.  Intuitively,  the 
aft  closure  flow  field  based  on  the  spherical  surface  flow  area  is  more  reas¬ 
onable  because  .e  velocity  vector  is  normal  to  the  flow  area;  this  is  not  the 
case  for  the  one  dimensional  flow  assumption.  The  convective  heat  transfer 
coefficients  for  both  flow  field  analyses  are  presented  in  Figure  3-9A. 

This  shows  that  the  modified  flow  analysis  resulted  in  lower  heat  transfer 
coefficients  in  the  aft  closure  only.  Two  dimensional  effects  in  the  blast 
tube  and  further  downstream  were  found  to  be  insignificant. 

Figure  3-11  illustrates  the  variation  of  the  ARGEIBL  heat  transfer 
coefficient  with  blast  tube  area  ratio  at  several  nozzle  locations.  Both  the 
5  and  the  18  percent  aluminum  propellants  exhibit  the  same  trends.  For  the 
blast  tube  entrance,  the  heat  transfer  coefficient  decreases  with  increasing 
blast  tube  area  ratio.  For  the  throat  and  the  exit  cone  entrance,  the 
convective  coefficient  increases  with  increasing  blast  tube  area  ratio. 

Figure  3-12  presents  the  heat  transfer  coefficient  distribution  for 
the  slotted  grain  port  for  a  nozzle  with  a  blast  tube  area  ratio  of  2.0  and 
for  a  5  percent  aluminum  propellant.  Triese  curves  apply  only  to  a  portion 
of  the  nozzle  which  is  in  the  same  plane  as  the  slot  in  the  propellant 
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calculate  the  free  stream  properties  throughout  the  nozzle  for  each  propellant 
These  properties  include  static  temperature,  static  pressure,  velocity,  density, 
modecular  weight  and  species  concentration  as  a  function  of  axial  distance  in  the 
nozzle.  This  program  also  calculates  some  of  the  more  significant  transport 
properties. 

In  calculating  the  free  stream  properties,  the  assumption  was  made  that 
local  thermodynamic  and  chemical  equilibrium  existed  throughout  the  nozzle  and 
that  the  propellant  gases  go  through  an  isentropic  expansion  from  the  chamber 
state.  This  allowed  a  series  of  closed  system  equilibrium  solutions  to  be  per¬ 
formed  at  the  various  axial  locations.  The  ACE  program  first  calculates  the 
entropy  at  the  chamber  state,  and  this  entropy  is  held  constant  in  subsequent 
solutions  at  different  locations.  The  major  results  of  this  analysis  are  pre¬ 
sented  in  Figure  3-6  for  the  5  percent  aluminum  CTPB  propellant  and  in  Figure  3-7 
for  the  18  percent  propellant.  These  results  are  based  on  a  chamber  pressure 
of  2800  psia. 

3.5  HEAT  TRANSFER  BOUNDARY  CONDITIONS 

After  the  boundary  layer  edge  or  free  stream  properties  were  calculated 
the  heat  transfer  boundary  conditions  were  determined.  These  boundary  conditions 
include  both  the  convective  and  radiative  heat  transfer  coefficients.  The  convec¬ 
tive  heat  transfer  coefficients  apply  during  a  motor  firing  (Section  3.5.1)  and  the 
radiative  heat  transfer  coefficients  apply  both  d  -.ring  the  firing  and  during  the 
motor  cooldown  (Section  3.5.2). 

3.5.1  Convective  Boundary  Conditions 

To  define  the  convective  boundary  conditions,  the  heat  transfer  coeff¬ 
icient  as  a  function  o.  boundary  layer  running  length  was  calculated.  This 
coefficient  was  defined  based  on  an  enthalpy  driving  potential  and  is  given  by 

q  =  p  u  Cu  (H  -  H  ) 

Mconv  e  e  H  R  ew 

where  1conv  “  Total  convective  heat  flux  to  the  material 

p.u  C„-  Convective  heat  transfer  coefficient  (  P  -  edge 

G  G  n  “ 

density,  ug  -  edge  velocity,  CH  -  Stanton  number) 

Hr  -  Recovery  enthalpy 

Hgw  -  Edge  enthalpy  at  the  wall  temperature. 

Aside  from  its  dependence  on  nozzle  location  and  chamber  pressure,  the  heat 
transfer  coefficient  varies  with  wall  temperature,  nozzle  geometry  (blast  tube 
area  ratio),  and  the  type  and  configuration  of  the  propellant  grain.  The  heat 
transfer  coefficients  dependence  on  these  quantities  is  presented  in  the 
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following  paragraphs  for  the  cylindrical  grain  and  keyhole  port  grains, 
respectively.  For  this  analysis,  a  chamber  pressure  of  2800  psia  was  used. 

The  heat  transfer  coefficient  distribution  in  nozzles  fired  with  a 
cylindrical  grain  propellant  were  calculated  accounting  for  the  effects  of 
1)  both  5  percent  and  1<3  percent  aluminum  propellant,  2)  blast  tube  area 
ratios  of  1.3  and  2.0,  3)  and  wall  temperatures  of  2700°  and  5400°B.  Figure 
3-8  illustrates  the  heat  transfer  coefficient  distribution  predicted  by  the  ARGEIBL 
code  for  the  5  percent  aluminum  propellant  for  blast  tube  area  ratios  of  1.3  and 
2.0.  In  each  figure  the  sensitivity  to  wall  temperature  is  shown.  Figure  3-9 
illustrates  the  heat  transfer  coefficient  for  the  18  percent  aluminum  propellant. 
Note  that  the  18  percent  aluminum  CTPB  yields,  in  general,  a  lower  heat  transfer 
coefficient  than  the  5  percent  propellant.  Also,  the  sensitivities  of  each  propel¬ 
lant  to  wall  temperature  are  similar.  For  the  case  of  the  1.3  blast  tube  area 
ratio,  an  analysis  was  performed  to  determine  the  sensitivity  of  the  heat  transfer 
coefficient  on  the  two  dimensional  flow  field  effects  in  the  aft  closure,  and  this 
analysis  is  described  in  the  following  paragraph. 

As  mentioned  previously  (Section  3.3),  a  one  dimensional  flow  field  analysis 
was  used  for  calculations  involving  the  cylindrical  port  propellant  cartridge.  To 
account  for  the  sharp  convergence  angle  of  the  aft  closure,  an  analysis  was 
performed  in  which  the  effective  flow  area  in  the  aft  closure  was  assumed  to 
be  a  spherical  surface  whose  center  is  the  apex  of  the  cone  created  by 
extending  the  aft  closure  surface  to  the  motor  centerline.  This  surface  is 
defined  by  the  schematic  drawing  presented  in  Figure  3-10.  Intuitively,  the 
aft  closure  flow  field  based  on  the  spherical  surface  flow  area  is  more  reas¬ 
onable  because  .e  velocity  vector  is  normal  to  the  flow  area;  this  is  not  the 
case  for  the  one  dimensional  flow  assumption.  The  convective  heat  transfer 
coefficients  for  both  flow  field  analyses  are  presented  in  Figure  3-9A. 

This  shows  that  the  modified  flow  analysis  resulted  in  lower  heat  transfer 
coefficients  in  the  aft  closure  only.  Two  dimensional  effects  in  the  blast 
tube  and  further  downstream  were  found  to  be  insignificant. 

Figure  3-11  illustrates  tha  variation  of  the  ARGEIBL  heat  transfer 
coefficient  with  blast  tube  area  ratio  at  several  nozzle  locations.  Both  the 
5  and  the  18  percent  aluminum  propellants  exhibit  the  same  trends.  For  the 
blast  tube  entrance,  the  heat  transfer  coefficient  decreases  with  increasing 
blast  tube  area  ratio.  For  the  throat  and  the  exit  cone  entrance,  the 
convective  coefficient  increases  with  increasing  blast  tube  area  ratio. 

Figure  3-12  presents  the  heat  transfer  coefficient  distribution  for 
the  slotted  grain  port  for  a  nozzle  with  a  blast  tube  area  ratio  of  2.0  and 
for  a  5  percent  aluminum  propellant.  Tnese  curves  apply  only  to  a  portion 
of  the  nozzle  which  is  in  the  same  plane  as  the  slot  in  the  propellant 


cartridge.  As  shown  by  Figure  3-12,  the  predicted  heating  is  increased 
substantially  in  the  aft  closure  by  the  keyhole  port  grain,  but  the  heating  in 
the  remainder  of  the  nozzle  is  reduced.  As  time  progresses  the  heat  transfer 
coefficient  is  decreased  in  the  aft  closure  and  increased  downstream.  The 
steady  state  value  shown  is  the  same  as  that  of  a  nozzle  with  a  cylindrical  grain 
propellant.  Figure  3-13  shows  the  heat  transfer  coefficient  both  in  the  slot 
plane  and  opposite  to  the  slot  plane  for  the  18  percent  Al-CTPB  propellant 
and  a  1.3  area  ratio  blast  tube.  Figure  3-14  normalizes  the  heat  transfer 
coefficient  predicted  for  the  5  percent  A1  propellant  keyhole  port  using  the 
circular  port  flow  for  the  aft  closure,  blast  tube  entrance,  throat  package, 
and  exit  cone  as  a  baseline.  These  results  were  used  to  determine  the  most 
severe  heating  conditions  to  be  studied  in  the  thermal  screening  analysis  which 
is  presented  in  the  following  section. 

As  mentioned  earlier  the  heat  transfer  coefficient  distributions  given 
in  the  figures  presented  in  this  section  were  calculated  by  the  ARGEIBL  code. 
Previous  experience  has  shown  that  the  heat  transfer  coefficient  predicted  by 
the  ARGEIBL  code  is  high  and  should  be  multiplied  by  0.75;  that  is 

peUeCH  =  0,75  {peueCH) ARGEIBL  (3-2) 


where  ^peUeCH^ ARGEIBL-  Convective  heat  transfer  coefficient  predicted  by 

the  ARGEIBL  Code. 


The  factor  of  0.75  has  been  used  by  Aerotherm  in  previous  thermal  analyses  of 
solid  rocket  motors  (References  3  and  4) .  For  the  transient  in-depth  thermal 
analysis  presented  in  subsequent  sections,  Equation  3-2  is  modified  further 
to  account  for  the  actual  chamber  pressure  history  by  using  the  relationship 

peUeCH  =  0,75  ^ 2800  ^peueCH^ ARGEIBL  (3-3: 


where  PQ  is  the  instantaneous  n.  red  (or  predicted)  chamber  pressure  in  psia 
units. 

3.5.2  Radiative  Boundary  Conditions 

To  model  the  radiative  boundary  conditions  during  nozzle  firings,  a 
parallel  plate  model  was  uted.  This  model  applies  to  aluminized 
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propellants  and  assumes  that  the  particle  laden  stream  of  combustion  products  is 
optically  thick  and  that  it  exchanges  radiant  energy  with  the  surface  as  if  the 
stream  and  wall  were  parallel  plates-  In  this  way,  multiple  reflections  between 
the  wall  and  stream  were  taken  into  account.  In  addition,  the  assumption  was  made 
that  both  the  stream  and  wall  behave  as  gray  bodies  and  that  they  emit  and  reflect 
radiant  energy  diffusely.  Based  on  the  above  assumptions,  the  net  radiant  heat 
flux  relation  is  given  as 

q  net  rad  =  Eeff(°T'‘  ~  oT^)  (3-4) 

where 

e  -  Effective  emmissivity  =  - i, -  (3-5) 

eff  -L  +  -L 1 1 

Ew  Ew 

ew  -  Wall  material  emissivity 

es  -  Particle  laden  stream  emissivity 

o  -  Stefan-Boltzmannn  constant 

Ts  -  Free  stream  (edge)  temperature 

Tw  -  Wall  temperature. 

To  determine  the  effective  emissivity  using  equation  (3-5) ,  the  stream 
emissivity  was  defined  as 


ts  =  1  -  exp  (-C^j-  pD) 


(3-6) 


where 

C 

n 

P 

D 


-  Empirical  constant  (0.808) 

-  Percentage  of  aluminum  loading 

-  Local  density  of  propellant  combustion  species 
(lb/ ft3) 

-  Local  beam  length,  usually  taken  as  the  diameter 

( in . ) . 


To  model  the  radiative  boundary  conditions  during  periods  of  nozzle  cooldown,  a 
different  model  was  used.  The  nozzle  surface  being  analyzed  was  considered  as 
exchanging  radiant  energy  with  only  one  source  which  was  the  environment  beyond 
the  exit  plane  of  the  nozzle.  This  implies  that  the  remaining  interior  nozzle 
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can  ba  approximated  as  an  isothermal  surface  which  is  at  the  same  temperature 
as  the  one  being  analyzed.  The  radiant  exchange  with  the  environment  beyond 
the  exit  plane  of  the  nozzle  was  based  on  the  assumption  of  diffuse  gray  body 
radiation.  Therefore,  the  surface  being  analyzed  would  radiate  energy  to  this 
environment  in  the  form  of 

q  net  rad  =  Few«j  (T^  -  T^)  (3-7) 

where 

F  -  View  factor  from  surface  to  environment  beyond  exit 

plane. 

Tg  -  Environment  (ambient)  temperature. 

If  T1"  is  neglected  in  relation  to  T"*  Equation  (3-7)  reduces  to 

cl  " 

q  net  rad  -  -Fe  oT11 
w  w 

3.6  MATERIAL  CHEMICAL  AND  THERMAL  PROPERTY  DATA 

This  section  summarizes  the  material  chemical  and  thermal  properties 
used  for  making  material  performance  predictions  using  the  CMA  and  ASTHMA 
computer  codes. 

The  thermal  properties  include 

•  density  -  virgin  material 

-  char  material 

•  material  specific  heat  -  virgin  material 

-  char  material 

•  material  thermal  conductivity  -  virgin  material 

-  char  material 

and  the  chemical  properties  include 

•  resin  mass  fraction 

•  effective  resin  and  reinforcement  molecules 

•  heat  of  formation  -  virgin  material 

-  char  material 

•  constants  in  Arrhenius  type  decomposition  note  equation  for  the 
resin  system 
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Material  thermal  properties  for  most  materials  were  obtained  from  the 
references  listed  in  Table  3-1.  Thermal  properties  for  R-155  were  not  avail¬ 
able,  thus  limiting  its  analysis  to  steady  state  calculations.  The  thermal 
properties  for  MX2625  were  estimated  (Figure  3-15)  from  those  of  MX2646 
Ref.  6).  The  thermal  conductivity  weighting  function  is  described  in  Table  3-2. 
The  thermal  properties  of  Pycobond  were  estimated  to  be  those  shown  in  Figure 
3-16  and  are  an  average  of  the  conductivities  of  ATJ  graphite  with  and  against 
the  grain.  Figure  3-17  gives  the  estimated  thermal  conductivity  of  Durez  16771-1 
which  was  obtained  by  decreasing  the  conductivity  of  MX2600  by  25%.  The  thermal 
conductivity  weighting  function  presented  in  Table  3-3  was  used.  The  nominal 
material  density,  resin  mass  fraction,  and  resin  residual  are  presented  in 
Table  3-4.  The  char  density  was  calculated  from  this  information.  Also 
presented  in  Table  3-4  are  the  chemical  formulas  for  the  resin  and  char  for 
each  material  except  R-155.  A  heat  of  formation  value  of  -1100  Btu/lb  was 
used  for  the  phenolic  resin  and  a  value  of  -6520  Btu/lb  was  used  for  the  SiOj. 

The  decomposition  kinetics  for  the  phenolic  resin  are  defined  in  Table  3-5. 

'tt'e  standard  chemical  information  was  not  available  in  the  literature 
for  R-155.  The  values  shown  in  Table  3-6  were  obtained  from  Reference  16, 
and  are  sufficient  for  calculating  both  the  steady  state  surface  thermochemistry 
map  and  the  steady  state  surface  recession  and  temperature. 
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TABLE  3-1 

SUMMARY  OF  MATERIAL  THERMAL 
PROPERTY  DATA 


Material 

Source  of 

Thermal  Properties 

MXCE  280 

Ref.  13 

MXSE  280 

Ref.  13 

MXSC  195 

Ref.  13 

MXi625 

Figure  3-15 

R-155 

M/A 

MX4926b 

Ref.  14 

FM5055 

Same  as  MX4926 

MX2600 

Ref.  14 

G-90 

Ref.  4 

KyroStranu 

Ref.  15 

Pycobond 

Figure  3-16 

Durez  16771-1 

Figure  3-17 

PG 

Ref.  17 

ATJ 

Ref.  18 

a  The  values  of  thermal  conductivity  used  for 
these  materials  are  presented  in  Section  5. 

b  Thermal  properties  suitable  for  the  ASTHMA 
code  are  presented  in  Section  5. 


TABLE  3-2 

MX2625  THERMAL  CONDUCTIVITY  WEIGHTING  FUNCTION 


Mass  Fraction 

Weighting  Function 
Virgin  Char 

1 .00 

1 .4 

.0 

.90 

.4 

.0 

.60 

.4 

.0 

.35 

.0 

1.0 

.0 

.0 

1.0 

TA8LE  3-3 

OUREZ  16771-1  THERMAL  CONDUCTIVITY  WEIGHTING  FUNCTION 
Mass  Function  Weighting  Function 


Virgin  Char 


1.00  1.0  .0 


.8  .6  .0 

.65  .5  .0 

.25  .0  1.0 

.0  .0  1.0 


TABLE  3-4 

MATERIAL  CHEMICAL  DATA  USED  FOR  THERMAL  SCREENING  ANALYSIS 
AND  IN-DEFTH  THERMAL  ANALYSIS 
AND  PRETEST  PREDICTIONS 


Material 

Nominal 

Density 

Lfc/Ft 

Nominal 

Resin 

Mass 

Fraction 

Resin 

Residual 

Assumed 

Resin 

Elemental 

Formula 

Peinforcement 

Elemental 

Formula 

MXCE  280,  Fibrous  carbon 
Impregnated  with  an  elas¬ 
tomeric  modified  phenolic 
resin 

78.0 

0.475 

0.40 

w 

C 

MXSE  280,  Fibrous  silica. 
Impregnated  with  an  elas¬ 
tomeric  modified  phenolic 
resin 

96.2 

0.350 

0.50 

w 

sio2 

MXSC  195,  Fibrous  poly¬ 
crystalline  carbon  silica 
Impregnated  with  phenolic 
resin 

96.8 

0.400 

0.45 

C6H6° 

C/Si02 

MX2625,  Fibrous  silica, 
Impregnated  with  phenolic 
resin 

111.8 

0.250 

0.5 

w 

sio2 

R-155,  ASB/Rubb r« 

79.5 

N/A 

N/A 

N/A 

N/A 

MX4926,  Carbon  cloth 
Impregnated  with  phenolic 
resin 

89.4 

0.345 

0.40 

w 

C 

MX2600 ,  Fibrous  silica, 
Impregnated  with  phenolic 
resin 

108.7 

0.315 

0.50 

w 

Si02 

G-90  Bulk  graphite 

117.5 

- 

- 

- 

C 

Pyrostrand,  graphite 
base  with  graphite  bond 

112.0 

- 

- 

- 

C 

Pycobond,  graphite 
fiber  with  grahpite  bond 

103.0 

0.40 

Durez  16771-1 

115.4 

0.345 

0.5 

c6h6o 

Si02 

FM5055,  Carbon 

Impregnated  with  phenolic 
resin 

89.4 

0.345 

0.40 

c6H6° 

C 

Pyrolytic  graphite 

137.3 

- 

C 

ATJ 

— 

108.0 

- 

- 

C 

igure  3-1.  Blast  Tube  Nozzle  Design 
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FIGURE  3-2.  LOADED  PROPELLANT  CARTRIDGE 
SLOTTED  GRAIN 


FIGURE  3-4.  VARIATION  OF  NON-DIMENSIONAL  FLOW  PARAMETER  WITH 

BOUNDARY  LAYER  RUNNING  LENGTH  FOR  THE  SLOTTED  GRAIN 
A)  BLAST  TUBE  AREA  RATIO  =  2.0,  SLOT  PLANE 


FORWARD 

ENTRANCE 


EFFECTIVE  FLOW  AREA  = 
^ - (WIDTH  OF  SLOT)  X  (r] 


VELOCITY  DIRECTION 


INSULATION 


FIGURE  3-5.  DEFINITION  OF  EFFECTIVE  FLOW  AREA  IN  THE 
AFT  CLOSURE  -  SLOTTED  GRAIN  CONFIGURATION 
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FIGURE  3-6.  (CONTINUED) 

B)  THROAT  REGION 
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FIGURE  3-8.  HEAT  TRANSFER  DISTRIBUTIONS  FOR  5  PERCENT 
AL-CTPB  PROPELLANT 

A)  BLAST  TUBE  CONTRACTION  RATIO  =1.3 


800  PSIA 


B)  BLAST  TUBE  CONTRACTION  RATIO  =  2.0 
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FIGURE  3-10.  AFT  CLOSURE  CROSS  SECTIONAL  AREA 
FOR  MODIFIED  FLOW  FIELD  ANALYSIS 
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FIGURE  3-11.  (CONCLUDFO) 
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FIGURE  3-12.  HEAT  TRANSFFR  DISTRIBUTIONS  FOR  THE 
SLOTTED  GRAIN  PORT  AT  VARIOUS  TIMES 
AFTER  START-UP 


FIGURE  3-13.  HEAT  TRANSFER  DISTRIBUTIONS  FOR  THE  SLOTTED  GRAIN 
PORT  AT  VARIOUS  TIMES  AFTER  STARTUP 

A)  IN  CIRCUMFERENTIAL  PLANE  WITH  SLOT 
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FIGURE  3-14.  EFFECT  OF  KEYHOLE  PORT  FLOW  ON  HEAT  TRANSFER  TO 
VARIOUS  LOCATIONS  OF  THE  BLAST  TUBE  NOZZLE  IPLANE 
OF  SLOT) 


FIGURE  3-15.  ESTIMATED  THERMAL  PROPERTIES  OF  MX2625  SILICA  PHENOLIC 
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FIGURE  3-16.  ESTIMATED  THERMAL  PROPERTIES  OF  PYCOBOND  410-65 
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THERMOCHEMICAL  SCREENING  ANALYSIS  AND  IN-DEPTH  THERMAL  ANALYSIS 

The  objective  of  the  thermochemical  screening  analysis  was  to  predict  the 
surface  recession  rate  and  surface  temperature  for  classes*  of  materials  when 
subjected  to  the  high  chamber  pressure  environment.  As  described  in  Section  2, 
this  screening  is  performed  by  coupling  of  the  ACE  program  to  a  steady  state 
energy  balance  solution  (SSA  computer  program) .  This  calculation  does  not 
require  a  knowledge  of  the  material  thermal  properties  because  of  the  steady 
state  assumption,  but  it  does  require  a  knowledge  of  the  material  chemical 
properties.  The  results  of  the  thermochemical  screening  analysis  provide  a 
basis  for  ranking  materials  and  also  provide  the  surface  boundary  conditions 
for  performing  a  conservative  in-depth  analysis.  The  parameters  considered  in 
the  analysis,  in  addition  to  material  class,  were  the  aluminum  loading  fraction, 
the  blast  tube  area  ratio,  the  propellant  grain  configuration,  and  the  nozzle 
location. 

In  performing  the  thermochemical  screening  analysis,  the  candidate  nozzle 
materials  being  considered  for  the  high  chamber  pressure  nozzle  were  divided 
into  categories  according  to  material  class.  In  addition,  a  typical  material 
was  selected  for  each  material  class.  The  selection  of  candidate  materials  and 
the  division  of  these  materials  into  material  classes  were  performed  jointly 
by  UTC  and  Aerotherm.  Also,  the  nozzle  location  for  each  material  class  and  the 
particular  material  within  each  class  which  was  considered  as  typical  were 
defined  by  this  joint  effort.  The  material  classes  considered,  the  nozzle 
location,  and  the  typical  material  within  each  material  class  are  presented  in 
Table  4-1.  In  Table  4-2  both  the  material  location  and  the  propellant  considered 
are  defined.  The  approach  taken  was  to  first  make  simplified  calculations  to 
determine  the  steady  state  surface  recession  and  surface  temperature  of  each 
material  to  both  the  5  and  18  percent  Al-CTPB  propellants  for  various  area 


* 

The  term  "material  class"  refers  to  materials  which  are  composed  of  similar 
reinforcements;  thus,  all  carbon  and  graphite  reinforced  ablation  materials 
would  form  one  class  and  all  silica  reinforced  materials  would  form  another 
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ratios  within  the  nozzle  section  of  interest.  Given  this  information,  an  in- 
depth  thermal  analysis  was  made  as  to  which  material,  propellant,  and  geometry 
provided  the  most  severe  thermal  response  in  each  of  the  three  nozzle  components 
of  interest:  t.  e  aft  closure,  the  blast  tube,  and  the  exit  cone.  These  analysis 
results  which  were  obtained  using  option  2  of  the  CMA  program  provided  a  basis 
for  establishing  the  material  thickness  requirements  in  the  nozzle  design. 

As  shown  in  Table  4-2,  a  transient  CMA  analysis  was  performed  for  the 
G-90  and  Pyrostrand  blast  tube  configurations.  A  transient  analysis  was  required 
for  these  materials  because  the  steady  state  assumption  does  not  apply  for  the 
following  reasons: 

»  The  surface  recession  is  kinetically  controlled. 

•  The  material  thermal  conductivity  is  relatively  high. 

Section  4.1  presents  the  surface  thermochemistry  data  and  the  heat  flux 
boundary  condition  information  obtained  from  Section  3  which  are  required  in 
order  to  perform  the  screening  and  in-depth  analysis.  Section  4.2  presents 
the  thermochemical  screening  analysis  results  and  Section  4.3  presents  more  de¬ 
tailed  in-depth  thermal  analysis.  Also  included  in  this  section  are  two  in-depth 
thermal  analyses  using  transient  convective  heating  (rather  than  steady  state) 
boundary  conditions  to  determine  the  thermal  and  ablative  response  of  the  G-90  and 
Pyrostrand  materials. 

4.1  SURFACE  THERMOCHEMISTRY  DATA  AND  HEA"  FLUX  BOUNDARY  CONDITIONS 

Figures  4-1  and  4-2  present  the  surface  thermochemical  responses  of  the 
charring  ablating  materials  considered  in  the  screening  analysis.  For  internally 
decomposing  materials  a  computer  code  such  as  CMA  calculates  material  response 
transients  by  solvinq  iteratively  a  surface  energy  balance  which  requires  a  para¬ 
metric  input  of  both  pyrolysis  off-gas  rates  and  char  removal  rates  as  a  function 
of  surface  temperature.  Surface  thermochemistry  information  of  this  type  is 
used  in  Section  5  for  transient  calculation  of  internally  decomposing  materials 
in  the  throat  package.  However,  for  preliminary  ranking  of  materials,  this 
degree  of  sophistication  is  ’.at  required,  and  thus  the  steady  state  assumption 
is  warranted.  For  the  steady  state  surface  response  of  a  charring  material, 
the  ratio  of  the  mass  removal  rates  of  the  pyrolysis  gas  and  char  satisfy  the 
steady  state  mass  conservation  equation.  Figures  4-1  and  4-2  are  thermochemistry 
maps  based  on  this  steady  state  assumption  in  which  the  non-dimensional  mass 
removal  is  a  combination  of  both  pyrolysis  gas  and  char  removal. 

For  the  cases  considered,  surface  response  maps  ai.e  functions  of  stream 
temperature,  pressure,  and  propellant.  However,  since  the  variation  of  the  maps  with 
temperature  and  pressure  is  small,  these  maps  were  calculated  for  only  one  state 


52 


point  for  each  material.  This  temperature  and  pressure  was  selected  as  the 
average  for  the  locations  at  which  each  material  was  considered.  As  a  result, 
two  thermochemistry  maps  were  calculated  for  each  material;  one  for  5  and  another 
for  18  percent  Al-CTPB  propellant.  These  maps  are  used  in  Section  4.2  to  per¬ 
form  the  thermochemical  screening  analysis. 

Figures  4-3  and  4-4  present  the  surface  thermochemistry  maps  for  the  two 
blast  tube  materials  for  which  in-depth  thermal  analyses  using  a  convective 
heating  boundary  condition  were  performed.  The  G-90  surface  response  was  cal¬ 
culated  for  a  blast  tube  pressure  of  160  atm  for  the  18  percent  Al-CTPB  pro¬ 
pellant.  The  pyrostrand  response  was  performed  for  a  pressure  of  162  atm  for 
both  the  5  and  18  percent  Al-CTPB  propellants. 

4.2  THERMOCHEMICAL  SCREENING  ANALYSIS 

The  thermochemical  screening  analysis  was  used  as  a  means  to  quickly 
evaluate  the  material  performance  of  a  large  number  of  materials  at  several 
nozzle  locations.  This  analysis  was  accomplished  by  solving  a  steady  state  sur¬ 
face  energy  balance  for  each  material,  geometry,  location,  and  propellant 
summarized  in  Table  4-3.  The  energy  jalanoe  calculations  were  performed  using 
the  SSA  computer  code.  Input  to  this  code  consists  basically  of  a  heat  trans¬ 
fer  coefficient,  edge  enthalpy,  and  surface  thermochemistry  tables  for  each 
solution.  These  quantities  were  obtained  from  Sections  3.5.1,  3.4,  and  4.1 
respectively.  Table  4-3  summarizes  the  results  of  these  calculations.  A 
schematic  representation  of  the  analysis  locations  is  presented  in  Figure  4-5. 

The  results  presented  in  this  table  are  discussed  in  more  detail  in  the  following 
paragraphs . 

As  shown  in  Figure  4-5,  analyses  weie  performed  for  the  aft  closure  for 
area  ratios  corresponding  to  9.0  and  22.0.  Figures  3-8  and  3-9  show  that  a 
contraction  ratio  of  9.0  (A  =  7.9  inches)  experiences  the  most  thermally  severe 
boundary  conditions  in  the  aft  closure.  Steady  state  solutions  for  four 
materials  for  both  5  and  18  percent  A1  propellants  are  given  in  Table  4-3A  for 
the  location  corresponding  to  an  area  ratio  of  9.0.  Further  studies  of  in- 
depth  thermal  response  analysis  at  this  location  was  not  practical,  however, 
due  to  the  close  proximity  of  the  forward  entrance  cap.  The  difficulty  in  this 
analysis  was  that  the  surface  recession  normal  to  the  surface  was  greater  than 
the  available  thickness  of  the  aft  closure  material,  i.e.,  the  local  receding 
surface  intersected  the  forward  entrance  cap  material.  For  this  reason  further 
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thermochemical  screening  was  performed  at  a  contraction  ratio  of  22.0,  and  the 
results  are  given  in  Table  4-3B.  The  location  at  a  contraction  ratio  of  22.0 
is  shown  in  Figure  4-5  to  be  the  section  at  which  the  aft  closure  material  is 
thickest.  As  shown  in  Table  4-23,  the  MXSE  280  material  exposed  to  the  5  per¬ 
cent  A1  propellant  exhibited  the  largest  recession  rate. 

A  further  stop  in  the  screening  analysis  was  made  by  calculating  the  effect 
of  the  keyhole  port  grain  propellant  on  the  aft  closure  surface  recession.  The 
results  presented  in  Table  4-3B  shows  that  this  grain  configuration  increased 
the  steady  state  recession  by  a  factor  of  five  at  conditions  representing  the 
initiation  of  firing  but  that  this  factor  diminished  to  nearly  the  steady 
state  condition  after  five  seconds  of  firing  as  illustrated  in  Figure  3-5  and 
Figure  3-12.  Based  on  these  »-esults,  the  thickness  of  the  aft  closure  material 
must  be  sized  or.  the  surface  recession  for  MXSE  280  material  and  for  a  slotted 
grain  5  percent  Al-CTPB  propellant.  For  this  reason,  an  in-depth  thermal 
analysis  of  this  configuration  is  presented  in  Section  4. 2. 2.1. 

Thermochemical  screening  analyses  were  performed  at  four  different  blast 
tube  area  ratios.  Area  ratios  of  1.3  and  2.0  were  specified  as  design  variables 
for  the  two  nozzle  configurations  of  interest  while  3.0  and  4.0  were  selected 
to  gain  information  concerning  the  effect  of  local  area  ratios  resulting  from 
accumulated  surface  recession  on  the  instantaneous  recession  rate.  Table  4-3C 
summarizes  the  results  of  this  analysis.  These  calculations  were  made  at  the 
blast  tube  entrance  because  Figures  3-8  and  3-9  show  that  this  location  ex¬ 
periences  the  maximum  heating  in  the  blast  tube. 

For  the  materials  considered  m  this  analysis,  variations  in  blast  tube 
area  ratio  affect  both  the  steady  state  surface  recession  and  surface  temperature. 
Figure  4-6  illustrates  the  fact  that  increasing  the  blast  tube  radius  will  de¬ 
crease  surface  recession  and  will  increase  or  hold  constant  surface  temperature. 
This  figure  also  shows  that  the  propellant  with  a  greater  aluminum  loading 
fraction  (18  percent)  decreases  surface  recession  and  increases  surface  tempera¬ 
ture  over  the  5  percent  Al-CTPB  pliant. 

The  choice  of  the  worst  case  combination  of  material,  propellant,  and 
blast  tube  area  ratio  for  the  in-depth  thermal  analysis  was  highly  qualitative. 
Table  4-3C  shows  that  for  area  ratios  of  either  1.3  or  2.0  the  silica  phenolic 
has  at  least  a  factor  of  2  greater  initial  surface  recession  rates.  He  -:ever  the 
greatest  thermal  penetration  is  obtained  for  a  material  which  has  a  relatively 
high  surface  recession  and  high  thermal  conductivity.  For  this  reason,  the  com¬ 
bination  of  the  MXSC  195  material  and  the  5  percent  Al-CTPB  propellant  were 
chosen  for  further  in-depth  thermal  analysis.  MXSC  195  has  the  relatively 
high  conductivity  of  a  cirbon  phenolic  while  exhibiting  the  relatively  greater 
recession  rates  of  a  silica  phenolic.  A  blast  tube  area  ratio  of  2.0  was  chosen 
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due  to  the  wall  thickness  being  less  than  for  a  1.3  area  ratio  blast  tube. 

In  screening  materials  for  the  exit  cone,  two  locations  (expansion  ratios 
of  2.1  and  6.0)  and  two  materials  (MX  4926  carbon  phenolic)  were  studied.  The 
results  are  summarized  in  Table  4-3D.  For  further  in-depth  thermal  analysis, 

MX  4926  for  an  expansion  ratio  of  6.0  and  for  the  5  percent  Al-CTPB  propellant 
was  chosen.  Again,  the  choice  was  highly  qualitative.  The  carbon  phenolic  has 
a  higher  conductivity  than  the  silica  phenolic  and  the  expansion  ratio  of  6.0 
had  a  thinner  section  of  material  than  that  at  an  expansion  ratio  of  2.1. 

4.3  IN-DEPTH  THERMAL  ANALYSIS 

The  thermochemical  screening  analysis  presented  in  the  previous  section 
provided  estimates  of  the  steady  state  surface  temperature  and  surface  recession 
rate  for  a  matrix  of  materials,  locations,  and  environments  in  the  high  pressure 
blast  tube  nozzle.  Both  surface  recession  and  in-depth  temperatures  are  impor¬ 
tant  parameters  for  design  considerations,  and  for  this  reason,  in-depth  calcula¬ 
tions  were  performed  and  are  presented  in  this  section.  The  analyses  were  made 
using  option  2  of  the  one-dimensiona1  CMA  computer  code. 

Section  3.6  presents  the  material  thermal  properties  used  in  the  in-depth 
thermal  analyses.  These  were  not  needed  in  the  previous  section  since  steady 
state  calculations  deal  solely  with  the  surface  phenomena.  Section  4.3.)  pre¬ 
sents  the  analyses  assuminq  steady  state  ablation  response.  Three  locations 
were  analyzed  (the  aft  closure,  blast  tube  and  exit  cone)  using  pseudo-transient 
surface  recession  and  temperature  boundary  conditions  calculated  from  steady 
state  quantities.  The  conditions  (material,  propellant  and  location)  were 
chosen  to  be  those  defined  in  Section  4.2  which  would  best  test  the  adequacy 
of  the  nozzle  design.  Section  4.3.2  presents  in-depth  thermal  analyses  of 
two  materials,  G-90  bulk  graphite  and  Pyrostrand  fibrous  PG,  which  were  con¬ 
sidered  for  use  in  the  blast  tube.  These  computations  were  performed  using 
option  1  of  the  CMA  computer  code. 

4.3.1  Analysis  Assuminq  Steady-State  Ablation  Surface  Response 

In  Section  4.2,  the  thermochemical  screening  analysis  was  presented. 

From  these  results,  a  material,  location  and  propellant  were  chosen  in  order 
to  define  the  adequacy  of  the  proposed  blast  tube  nozzle  design  in  each  of  the 
three  nozzle  sections;  aft  closure,  blast  tube  and  exit  cone.  These  cases  are 
summarized  as  follows: 

•  Aft  closure  -  Section  4. 3. 2.1 

MXSE  280  silica  phenolic 
8  percent  Al-CTPB  propellant 
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22.0  contraction  area  ratio 
slotted  grain  configuration 

•  Blast  Tube  -  Section  4. 3. 2. 2 

MXSC  195  carbon/silica  phenolic 
5  percent  Al-CTPB  propellant 
2.0  contraction  area  ratio 

•  Exit  Cone  -  Section  4. 3. 2. 3 

MX  4926  carbon  phenolic 
5  percent  Al-CTPB  propellant 
6.0  expansion  area  ratio 

The  separate  analysis  for  each  location  is  presented  in  the  following 
three  sections. 

4. 3. 1.1  Aft  Closure  Analysis 

The  in-depth  analysis  in  the  aft  closure  was  performed  using  the  one- 
dimensional  CMA  code  which  requires  as  input  the  geometry,  the  time  dependent 
boundary  conditions  and  material  chemical  and  thermal  data.  The  one-dimensional 
in-depth  geometry  at  a  contraction  ratio  of  22.0  consists  of  1.32  inches  of 
MXSE  280  and  0.43  inches  of  steel.  The  material  property  information  was 
obtained  from  Section  3.6.  The  time  dependent  boundary  conditions  ware  derived 
as  mentioned  previously,  from  the  steady  state  solutions  giver,  in  Table  4-2B. 

Figure  4-7A  shows  the  estimated  pseudo-transient  surface  recession  bound¬ 
ary  condition  used  in  the  in-depth  calculation.  Figure  4-7B  shows  the  surface 
temperature  history  which  was  used  as  a  boundary  condition.  This  was  also  taken 
from  the  steady  state  results  presented  in  Table  4-3B. 

The  results  of  the  aft  closure  analysis  are  presented  in  Figures  4-8  and 
4-9.  The  former  shows  the  computed  in-depth  thermal  profiles  during  the  exposure 
and  during  cooldown.  The  profiles  during  the  motor  firing  which  are  shown  in 
Figure  4-8A  illustrate  the  high  surface  recession  coupled  with  the  relatively 
low  thermal  conductivity  common  silica  phenolics.  This  results  in  a  very  small 
thermal  penetration  (about  0.10  inch)  as  the  material  recedes.  The  in-depth 
profiles  calculated  for  the  cooldown  are  shown  in  Figure  4-8B  and  indicate  that 
the  thermal  penetration  depth  is  only  about  one  half  of  the  remaining  MXSE  280 
thickness  after  60  seconds  of  elapsed  time.  The  steel  casing  was  predicted  to 
have  a  negligible  temperature  rise  as  shown  in  Figure  4-8B.  Figure  4-9 
illustrates  the  surface  and  char  line  recessions  in  relation  to  the  total  material 
thickness. 
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4. 3. 1.2  Blast  Tube  Analysis 

The  blast  tube  analysis  was  performed  at  a  contraction  ratio  of  2.0  and 
a  location  which  corresponded  to  the  entrance  of  the  blast  tube.  The  one¬ 
dimensional  geometry  consisted  of  0.80  inch  of  MXSC  195,  0.79  inch  of  silica 
phenolic  backup  and  0.32  inch  steel.  The  material  property  information  was 
obtained  from  Section  3.6. 

The  time  dependent  boundary  conditions  were  calculated  from  the  steady 
state  surface  recession  and  surface  temperature  given  in  Table  4-2C.  These 
values  were  72.0  mils/second  and  3600°R,  respectively.  Figure  4-10A  shows  the 
calculated  pseudo-transient  surface  recession  history  used  for  the  in-depth 
thermal  analysis.  This  was  derived  by  assuming  the  blast  tube  recession  in¬ 
itially  was  at  its  steady  state  value  of  72.0  mils/second.  As  the  surface  recedes 
with  time  the  internal  diameter  gives  a  second  point  on  the  pseudo-transient 
curve.  This  method  was  continued  m  this  time  step  manner  to  give  tne  recession 
history  shown.  The  surface  temperature  boundary  condition  as  shown  in  Figure 
4-10B  remains  constant  at  the  steady  state  value  of  3600°R. 

Applying  the  above  information  to  the  1-D  CMA  analysis  led  to  the  results 
show.i  in  Fiqures  4-11  and  4-12.  Figure  4-11A  depicts  the  in-depth  thermal  re- 
sp_>nse  of  MXSC  195  to  5  percent  Al-CTPB.  The  relatively  high  conductivity  of 
tne  carbon/silica  phenolic  would  yield  a  large  thermal  penetration  depth  for  a 
slowly  receding  surface.  However,  the  MXSC  195  also  exhibited  such  a  severe 
surface  erosion,  the  thermal  penetration  was  kept  to  a  minimum  (about  0.10  inch). 
Although  local  thermal  penetration  was  small,  79  percent  of  the  material  was 
predicted  to  ablate.  Figure  4-1113  shows  that  the  larye  silica  p,  enolic  back¬ 
up  was  adequate  for  insulating  the  steel  structure.  The  predicted  surface  and 
char  line  recessions  are  shown  m  Figure  4-12. 

4. 3.1. 3  Exit  Cone  Analysis 

The  exit  cone  analysis  was  performed  at  an  expansion  area  ratio  of  6.0 
which  corresponds  to  the  exit  plane  of  the  exit  cone.  The  one-dimensional 
qeometry  consisted  of  an  0.72  inch  layer  of  MX  4926  carbon  phenolic  and  an  0.09 
inch  thick  backing  of  steel.  The  material  property  information  used  as  input 
into  the  CMA  code  was  given  in  Section  3.6. 

The  time  dependent  boundary  conditions  were  based  on  the  values  given 
in  Table  4-3D.  The  steady  state  surface  recession  is  given  as  10.2  mils/second 
and  the  surface  temperature  as  2738°R.  Similar  to  the  blast  tube  calculation, 
the  exit  cone  boundary  conditions  were  derived  assuming  the  transient  response 
consisted  of  a  series  of  steady  state  conditions  spaced  at  finite  time  steps 
during  the  firing.  This  surface  rate  history  is  given  in  Figure  4-13A  and  the 
surface  temperature  is  shown  in  Figure  4-1 3B. 
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Figure  4-14A  illustrates  the  in-depth  temperature  profiles  obtained  using 
the  CMA  code.  Due  to  the  relatively  mild  environment  at  the  exit  plane  of  the 
exit  cone,  the  surface  recession  was  low  compared  to  the  original  thickness  of 
the  MX  4926.  The  recession  was  estimated  to  be  17  percent  of  the  combined 
material  thicknesses.  The  thermal  penetration  beneath  the  receding  surface 
was  greater  in  this  case  than  in  both  the  previous  cases.  This  was  due  to  the 
higher  conductivity  of  the  carbon  phenolic  as  compared  to  the  matei lals  considered 
upstream  of  the  throat,  figure  4-llB  shows  that  in-depth  heating  is  felt  at  the 
steel  casing  at  about  30  seconds.  After  60  seconds  of  elapsed  time,  the  steel 
shell  is  predicted  to  have  a  temperature  rise  of  about  65°R.  Figure  4-12  shows 
the  predicted  surface  and  char  recessions  during  tne  firing. 

4.3.2  Analysis  Using  Convective  Heating  Boundary  Condition  to 

Study  Blast  Tube  Thermal  and  Ablative  Response 

This  section  presents  the  additional  one-dimensional  analysis  performed 
for  two  materials  in  the  blast  tube.  The  CMA  code  was  used  with  a  convective 
heating  boundary  condition  (option  1).  Variables  considered  are  as  follows: 

•  G-90  -  Section  4. 3. 3.1 

18  percent  Al-CTPB 

contraction  area  ratio  =  1.3 

blast  tube  entrance  and  exit  locations 

•  Pryostrand  -  Section  4. 3. 3. 2 

5  percent  and  18  percent  Al-CTPB 
contraction  area  ratio  -  1.3 
blast  tube  entrance  location 

Section  4. 3. 2.1  discusses  the  results  obtained  from  analyzinn  G-90  bulk  graphite 
and  Section  4. 3. 2. 2  presents  the  results  for  analyzing  a  Pyrostrand  blast  tube 
configuration . 


4. 3.2.1  G-90  Blast  Tube  Analysis 

The  CMA  code  was  used  analyze  the  in-depth  thermal  response  of  G-90 
at  two  blast  tube  locations.  The  one-dimensional  geometry  consisted  of  a  0.45 
inch  thick  G-90  sleeve  backed  by  0.45  inches  of  silica  phenolic  (cloth  re¬ 
inforcement  in  0°  layup  orientation)  and  0.075  of  steel  casing.  The  chemical 
and  thermal  properties  for  G-90  were  obtained  from  Section  3.6.  The  surface 
thermochemical  response  of  G-90  to  18  percent  Al-CTPB  propellant  was  calculated 
using  the  ACE  computer  code  and  is  presented  in  Figure  4-3.  This  thermochemistry 
map  was  used  for  both  the  entrance  and  exit  locations  of  the  blast  tube. 


Since  the  analysis  uses  a  convective  heating  boundary  condition,  the 
recovery  enthalpy,  incident  radiant  heat  flux  and  convective  heat  transfer  co¬ 
efficient  as  a  function  of  time  were  determined.  For  both  locations  in  the 
blast  tube,  the  recovery  enthalpy  and  incident  radiant  flux  were  calculated  to 
be  942.5  Btu/lb  and  865.0  Btu/ftJ-second ,  respectively.  The  time  dependent 
heat  transfer  coefticient  was  calculated  from  the  information  in  Figure  3-9A. 

This  figure  shows  that  for  a  constant  chamber  pressure,  the  predicted  heat 
transfer  coefficient  at  the  entrance  of  the  blast  tube  (>  ■  0.94  ft)  is  2.77 

it/ftJ-S"'~  and  at  the  exit  (*  -  1.54  ft)  is  lb/ft 'sec  and  that  it  decreases  along 
the  length  of  the  blast  tube.  These  coefficients  were  estimated  for  a  wall  temp¬ 
erature  of  5000°R  and  they  were  applied  to  Eqs.  3-2  and  3-3  to  obtain  the  heat 
transfer  coefficients  shown  in  Figure  4-16. 

Figure  4-17  shows  the  m-depth  thermal  profiles  for  the  G-90  blast  tube 
entrance.  After  11  seconds  of  exposure,  the  thermal  penetration  depth  was  pre¬ 
dicted  to  be  at  about  the  midpoint  of  the  silica  phenolic  backup.  After  60 
seconds,  a  temperature  rise  of  60°R  was  felt  at  the  steel  casing.  Figure  4-18 
illustrates  the  corresponding  thermal  profiles  found  at  the  exit  of  the  blast 
tube.  Conditions  at  this  location  were  less  severe  (lower  heat  transfer  coeffi¬ 
cient)  tnan  those  at  the  entrance.  For  this  reason  in-depth  temperatures  and 
thermal  penetrations  weie  slightly  less  at  the  exit  than  at  the  entrance. 

However,  for  both  locations,  the  temperature  rise  of  steel  structural  member  was 
approximately  the  same.  At  the  exit  of  the  blast  tube  the  steel  housing  was 
predicted  to  rise  58°R  in  temperature;  2°R  less  than  the  housing  at  the  en¬ 
trance.  The  surface  recession  is  presented  in  Figure  4-19A  for  both  cases. 

The  entrance  was  calculated  to  recede  0.24o  inch  and  the  exit  0.194  inch. 

Figure  4-19B  compares  various  temperature  histories  for  the  two  locations. 

4. 3. 2. 2  Pyrostrand  Blast  Tube  Analyses 

The  one-dimensional  CMA  code  was  used  to  analyze  the  in-depth  thermal 
response  of  a  Pyrostrand  sleeve  at  the  blast  tube  entrance  with  a  backup  of 
castable  carbon.  Calculations  were  made  for  both  5  and  18  percent  Al-CTPB 
propellants  for  a  blast  tube  contraction  ratio  of  1.3.  The  one-dimensional 
geometry  consisted  of  a  0.442  inch  thick  Pyrostrand  sleeve  backed  by  0.108 
inches  of  castable  carbon.  Chemical  and  thermal  properties  used  for  Pyrostrand 
are  presented  in  Section  3.6.  The  kinetically  controlled  surface  response  of 
Pyrostrand  was  illustrated  in  Figure  4-4.  Reference  2  .as  used  as  a  source  of 
the  surface  kinetics  constants,  These  constants  were  those  reported  for  layer 
oriented  pyrolytic  graphite. 
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Other  pertinent  calculations  which  were  input  to  the  transient  in-depth 
code  are  summarized  in  Table  4-4.  The  time  dependent  heat  coefficients  were 
calculated  using  Eqs.  3-2  and  3-3  and  are  presented  in  Figure  4-20.  This 
illustration  shows  that  the  heat  transfer  coefficient  at  the  entrance  of  the 
blast  tube  is  higher  than  that  at  the  exit. 

in-depth  thermal  profiles  due  to  exposure  to  both  5  percent  and  18  per¬ 
cent  Al-CTPB  are  shown  in  Figure  4-21.  In  general,  18  percent  Al-CTPB  results 
in  higher  surface  temperatures,  a  greater  surface  recession,  and  greater  thermal 
penetration  than  5  percent  Al-CTPB.  After  an  elapsed  time  of  60  seconds,  the 
Pyrostrand/castable  carbon  interface  was  predicted  to  be  at  a  temperature  of 
1 440°R  for  the  5  percent  Al-CTPB  propellant.  Using  18  percent  Al-CTPB,  the  inter¬ 
face  is  shown  to  be  at  a  temperature  of  1600°R.  Figure  4-22A  compares  the  surface 
recessions  for  both  cases.  The  18  percent  Al-CTPB  propellant  results  in  a  fac¬ 
tor  of  5  increase  in  surface  recession  as  compared  to  the  5  percent  Al-CTPB  en¬ 
vironment.  However,  the  total  surface  recession  amounts  to  less  than  two  per¬ 
cent  of  the  original  pyrostrand  thickness  in  both  cases.  These  low  surface 
recessions  are  characteristic  of  pyrolytic  graphite  materials  oriented  in  a 
layer  configuration.  Figure  4-22B  presents  the  surface  and  interface  tempera¬ 
ture  calculated  for  both  cases. 
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TABLE  4-1 


NOZZLE  LOCATION,  MATERIAL  CLASSES,  AND 
TYPICAL  MATERIALS  FOR  THE  HIGH  CHAMBER 
PRESSURE  SOLID  ROCKET  MOTOR 


|  Nozzle 
|  Location 

Material 

Class 

Typical 

Material 

Aft  Closure 

Silica  Phenolic 

MXSE  280 

Elastomeric 

R-155 

Carbon  Phenolic 

MXCE  280 

Carbon/Silica  Phenolic 

MXSC  195 

Blast  Tube 

Silica  Phenolic 

MX  2625 

MXSE  280 

Carbon  Phenolic 

MXCE  280 

Carbon/Silica  Phenolic 

MXSC  195 

Bulk  Graphite 

G-90 

Pyrolytic  Graphite  Fiber 

Pyrostrand 

Exit  Cone 

Silica  Phenolic 

MX  2600 

Carbon  Phenolic 

MX  4926 
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TABLE  4-2 


Pyrosti 


5  refers  to  5  percent  AL-CTPB  and  18  refers  to  18  percent  AL-CTPB 
Transient  CMA  option  1  analysis  performed  for  thes*  materials. 
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WALL  TEMPERATURE 


DIMENSIONAL  SURFACE  REMOVAL 


Figure  4-3.  Kinetically  Controlled  Surface  Response  of  G-90 
Bulk  Graphite  at  Blast  Tube  (e  =  1.3)  Conditions 
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Figure  4-4.  Kinetically  Controlled  Surface  Response  of  Pyrostrand 
Layer  Graphite  at  Blast  Tube  (c  “  1.3)  Conditions 


qure  4-5.  Blast  Tube  Nozzle  Design 
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Fiqure  4-11.  Blast  Tube  Temperature  Profiles 

A)  Specified  Wall  Temperature  and 
Recession  P.ate 
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figure  4-13.  Exit  Plane  In-Depth  Analysis  Boundary  Conditions 
A)  Recession  Rate  History 
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Figure  4—13.  Concluded. 

B)  Surface  and  Computed  Backwall 
Temperature  Histories 
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Figure  <-14.  Concluded 

B)  During  Cooldown 
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Figure  H..  ^“-"^Transfe,  Coefficient  Versus  Time 
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Temperature  Profiles  for  C-90  at 
Blast  Tube  entrance 

A)  During  Exposure 
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Figure  4-20. 


Convective  Heat  Transfer  Coefficient  Versus  Time  for 
Pyrostrand  Blast  Tube 
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Figure  4-21.  Concluded 

B)  Exposed  to  18  Percent  AL-CTPB  Propellant 
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Figure  4-22.  Response  Histories  for  Pyrostrand  Blast  Tube 


A)  Surface  Recession 
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SECTION  5 


NOZZLE  THROAT  PACKAGE  THERMAL  RESPONSE  PREDICTIONS 

The  tbarmochemieal  screening  and  in-depth  thermal  analyses  presented 
in  the  prev/.ous  sections  provided  a  rapid  but  approximate  estimate  of 
the  thermochemical  performance  of  the  candidate  materials.  To  obtain  more 
accurate  predictions,  in-depth  one-dimensional  analyses  were  performer'  through¬ 
out  the  nozzle  using  convective  heating  boundary  conditions.  Tnis  method  is 
considerably  more  accurate  in  that  it  couoles  the  materials  response  to  its  en¬ 
vironment  for  transient  conditions.  These  one-dimensional  in-depth  analyses  are 
suitable  for  modeling  most  materials  throughout  a  typical  nozzle  where  two- 
dimensional  conduction  effects  are  negligible.  One  region  in  which  two-dimensional 
effects  are  important  is  the  throat  package.  In  this  section,  the  throat  package 
of  the  nozzle  with  a  blast  tube  area  ratio  of  1.3  and  an  18  peicent  Al-CTPB 
propellant  is  analyzed.  The  ASTHMA  computer  code  was  used  to  assess  the  two- 
dimensional  effects,  and  the  CMA  (one-dimensional)  code  was  used  at  different 
locations  for  comparison. 

The  input  to  the  ASTHMA  code  is  similar  to  that  used  by  the  CMA  code. 
Propellant  data,  the  flow  field  analysis,  free  stream  properties,  and  the  heat 
transfer  boundary  conditions  were  all  obtained  from  the  information  presented  in 
Section  3.  Section  5.1  presents  the  throat  package  geometry  and  Section  5.2 
summarizes  the  surface  thermochemistry  and  material  property  data.  The  thermal 
response  predictions  for  the  throat  package  are  given  in  Section  5.3. 

5.1  THROAT  PACKAGE  GEOMETRY 

Figure  5-1  illustrates  the  throat  package  configuration  which  was  analyzed. 

It  consists  of  six  pyrolytic  graphite  (PG)  washers  backed  by  a  section  of  ATJ 
graphite  which  in  turn  is  backed  by  silica  phenolic.  The  entrance  cap  is  made 
of  G-90  bulk  graphite  and  the  retaining  ring  consists  of  MX4926  carbon  cloth 
phenolic.  Two  teflon  rings  were  inserted  between  three  of  the  PG  washers  to 
allow  for  thermal  expansion. 

5.2  SURFACE  THERMOCHEMISTRY  AND  MATERIAL  PROPERTY  DATA 

As  mentioned  in  Section  2,  the  heat  conduction  solution  for  a  nozzle  material 
using  the  ASTHMA  or  CMA  codes  require  a  parametric  input  of  surface  removal  rate 
as  a  function  of  wall  temperature  and  pressure  for  non-decomposing  materials. 
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These  calculat  tons  wer.  done  using  the  ACE  computer  code  and  are  presented  in 
Figure  5-2  fo;  all  the  materials  in  the  throat  package.  Note  that  the  PG 
kinetic  curves  are  based  on  "old"  kinetics.  Differences  resulting  from  the 
"old  •  kinetics  are  presented  in  Section  5.3.2.  For  materials  which  exhibit 
internal  decomposition,  the  surface  mass  removal  rate  is  also  dependent  on  the 
pyrolysis  off  gas  rate.  Since  the  ASTHMA  code  does  not  model  internal  decompo¬ 
sition  (Se-  ticn  2),  a  thermochemistry  solution  assuming  steady  state  ablation 
was  used  for  the  carbon  phenolic  retaining  ring,  and  this  assumption  eliminates 
the  dependence  of  surface  removal  rate  on  the  pyrolysis  off-gas  rate.  This 
solution  is  also  shown  in  Figure  5-2.  In  a  latter  section  (Section  5.3.2)  it 
will  ba  shown  that  this  steady  state  assumption  introduces  little  error  into 
the  total  surface  recession  calculations  when  compared  to  a  one-dimensional 
CMA  solution.  The  surface  thermochemistry  calculations  for  the  transient  one- 
dimensional  solution  are  shown  in  Figure  5-3. 

For  edge  PG  and  the  two  bulk  graphites,  the  source  of  the  thermal  proper¬ 
ties  were  defined  in  Section  3.  The  thermal  properties  for  edge  PG  and  ATJ  are 
presented  in  Figures  5-4  and  5-5.  Since  ASTHMA  doer  not  account  for  internal 
decomposition,  the  thermal  properties  for  the  carbon  and  silica  phenolics  were 
modified  to  accomodate  this.  The  thermal  conductivities  are  weighted  averages 
between  those  of  virgin  and  char  materials  and  the  specific  heats  account  for 
energy  absorbed  by  the  material  diring  decomposition.  For  this  reason  the 
thermal  properties  for  MX  4926  carbon  phenolic  and  the  FM5504  silica  phenolic 
are  different  for  the  CMA  analysis  than  for  the  ASTHMA  2-D  analysis.  These 
properties  were  obtained  from  Reference  19  and  are  presented  in  Figures  5-6  and 
5-7.  The  chamber  pressure  histories  used  for  the  analysis  were  simplified  by 
UTC  and  are  presented  in  Figure  5-8. 

5.3  V1ERMAL  RESPONSE  PREDICTIONS 

"sing  (^e  information  presented  in  the  previous  sections,  predictions 
were  made  to  er.t  imate  the  thermal  response  of  the  components  of  the  throat 
package.  A  two-dimensional  analysis  performed  using  the  ASTHMA  code  is  presented 
in  Section  5.3.1.  One-dimensional  CMA  analyses  were  done  for  the  throat  and 
MX  4926  retaining  ring.  These  results  and  their  comparisons  with  the  2-D  analysis 
results  are  given  in  Section  5.3.2.  Section  5.3.3  presents  an  analysis  of  the 
throat  section  assuming  an  insulated  boundary  condition  at  the  backwall  of  the 
PG  washers. 
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5.3.1  Two-Dimensional  Analysis 

Using  the  geometry,  thermochemistry  solutions,  and  material  properties 
given  previously,  the  ASTHMA  code  was  used  to  predict  surface  temperatures,  in- 
depth  temperatures  and  surface  recession  for  an  11  second  exposure  to  18  percent 
A2-CTPB  propellant. 

Figure  5-9  presents  the  predicted  isothermal  lines  generated  within  the 
throat  package  for  several  elapsed  times  during  the  exposure  and  cooldown. 

These  plots  illustrate  how  the  high  radial  conductivity  of  the  PG  washers  allows 
energy  to  dissipate  from  the  throat  surface  to  the  backup  materials.  Also 
shown  are  the  insulating  effects  of  the  low  conductivity  carbon  phenolic  re¬ 
taining  ring.  Predicted  surface  recession  is  also  indicated.  In  this  case  the 
PG  washers  were  predicted  to  recede  6.3  mils  at  the  throat,  an  average  of  7.7 
mils  upstream  of  the  throat,  and  an  averaqe  of  2.5  mils  downstream  of  the 
throat.  The  MX4926  exhibited  the  most  erosion;  ranging  from  .225  inches  near 
the  PG  washers  to  .165  inches  at  the  downstream  end  of  the  insert.  The  G-90 
bulk  graphite  was  predicted  to  recede  approximately  .160  inches. 

In  calculating  the  heat  transfer  coefficient  distribution  use  for  this 
analysis  the  probable  separation  of  the  turbulent  boundary  layer  which  occurs 
downstream  of  the  PG  washers  vas  not  accounted  for.  It  was  assessed  that 
a  separation  of  the  flow  and  suosequent  reattachment  on  the  MX4926  insert  would 
result  in  a  higher  than  calculated  heat  transfer  rate.  Figure  5-10  shows  the 
surface  contour  and  isothermal  lines  in  the  throat  region  after  the  exposure 
period  with  an  estimate  of  the  gouging  caused  by  the  increased  heat  transfer. 

All  of  the  isotherm  figures  illustrate  the  estimated  erosion  of  the  PG  washer 
caused  when  the  upstream  G-90  insert  recedes  faster  than  the  washer.  This 
estimate  is  necessary  since  the  ASTHMA  code  accounts  only  for  environmental 
heating  on  surfaces  parallel  to  the  nozzle  centerline. 

5.3.2  One-Dimensional  Analysis  Results  and  Their  Comparison  to  the  Two- 

Dimensional  Analysis 

Two  locations  in  the  throat  package  were  chosen  to  perform  a  one-dimensional 
analysis  for  comparison  to  the  2-D  results.  These  locations  are  at  the  throat 
and  at  the  expansion  ratio  of  1.61  which  is  in  the  region  of  the  MX  4926  throat 
retaining  ring. 
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A  total  of  four  one-dimensional  analyses  were  performed  at  the  tnroat. 

These  correspond  to 

•  18  percent  A1-CTP3,  "old"  kinetics 

m  5  percent  Al-CTPB,  "old"  kinetics 

«  18  percent  Al-CTPB,  "new”  kinetics 

5  percent  Al-CTPB,  "new"  kinetics 

The  "new"  kinetics  correspond  to  the  ones  currently  in  the  GASKET  code  and 
documented  m  Rererence  2. 

Figure  5-11  shows  the  in-depth  temperature  histories  for  both  the  5  percent 
and  18  percent  Al-CTPB  propellant  cases  using  the  1-D  analysis  and  the  18  percent 
Al-CTPB  case  using  the  2-D  analysis.  The  thermal  profiles  a:.t*  shown  at  the 
conclusion  of  the  motor  firing  and  at  two  times  during  cooldown  and  are  for  the 
"old"  kinetics.  These  results  show  that  the  in-depth  thermal  response  given  by 
the  one  and  two  dimensional  analyses  are  almost  identical.  This  is  due  to  the 
highly  anisotropic  thermal  conductivity  <^f  the  PG  washers  which  allows  relatively 
little  energy  to  flow  in  the  axial  direction.  Figure  5-11  also  illustrates  the 
efficiency  of  the  silica  phenolic  backup  ns  an  insulator.  After  60  seconds  from 
startup,  the  m  alysis  shows  that  the  steei  shell  is  essentially  at  room  tempera¬ 
ture.  This  figuia  also  shows  that  the  18  percent  Al-CTPB  propellant  results  i  i 
considerably  higher  in-depth  and  surface  temperatures  than  the  5  percent  Al-CTPB 
propellant.  This  is  &  direct  result  of  the  higher  flame  temperature  of  the 
Al-CTPB  propellant. 

The  corresponding  1-D  analysis  for  both  of  the  propellants  using  the  "new" 
kinetics  showed  little  effect  in  the  surface  and  in-depth  temperature  profiles. 

For  the  18  Al-CTPB  propellant,  the  "new"  kinetics  led  to  only  a  50°R  maximum 
drop  in  surface  temperature  during  the  firing.  The  5  Al-CTPB  propellant  with  the 
"new"  kinetics  resulted  in  only  20°R  maximum  drop  in  temperature.  Histories 
for  all  four  1-0  cases  and  the  2-D  case  are  presented  in  Figure  5-12.  Although 
the  "new"  kinetics  lead  to  similar  temperatures,  the  difference  in  predicted 
surface  recession  is  significant. 

The  predicted  surface  recessions  for  the  throat  are  shown  in  Figure  5-13. 

In  general  the  predicted  surface  erosion  was  on  the  order  of  3  to  4  times  greater 
for  the  "new"  kinetics  than  for  the  "old"  kinetics. 
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An  expansion  ratio  of  1.61  was  chosen  to  perform  one-dimensional  CMA  analyses 
for  the  MX4926  throat  retaining  ring.  This  location  is  shown  in  Figure  5-9A. 

These  results  were  used  to  compare  1-1?  2-D  results  and  to  determine  the  accuracy 
of  the  ASTHMA  code  in  handling  charring  ablators.  Figure  5-1*  presents  the 
in-depth  temperature  profiles  for  both  the  one  and  two  dimensional  solutions  after 
an  11  second  firing.  The  CMA  analysis  gave  lower  temperatures  and  a  lower  surface 
recession.  Since  the  CMA  code  was  designed  to  analyze  internally  decomposing 
materials,  it  was  felt  that  the  surface  recession  and  temperatures  close  to  the 
surface  were  predicted  more  accurately  by  this  code.  Actual  2-D  effects  near 
the  surface  were  considered  negligible  at  the  analysis  location  as  evidenced  by 
the  isothermal  lines  shown  in  Figure  5-9.  This  figure  also  indicates  that  2-D 
heat  conduction  effects  are  sizable  at  deptns  greater  than  the  PG  washer  backwall 
for  an  elapsed  time  of  11  seconds  at  tl.‘.  location  of  interest.  Physically,  the 
IG  washers  conduct  heat  into  the  ATJ  rraph-te  backup  which  causes  axial  he< t 
conduction  into  the  neighboring  MX4926  retainer.  For  these  reasons,  the  actual 
in-depth  thermal  profile  is  a  combination  of  the  profiles  given  by  the  two  analyses. 
This  estimate  is  shown  in  Figure  5-14. 

In-depth  thermal  profiles  predicted  by  the  2-D  analysis  for  six  different 
times  during  the  motor  firing  and  subsequent  cooldown  are  illustrated  in  Figure  5-15. 
Th*.i<  figure  indicates  that  the  thermal  protection  provided  by  the  silica  phenolic 
backup  is  sufficient  to  keep  the  steel  shell  from  a  significant  rise  in  temperature, 
"igures  5-16  and  5-17  compare  the  surface  recessions  and  surface  temperatures 
given  by  the  1-D  and  2-D  analyses.  As  mentioned  above,  the  one-dimensional 
analysis  gives  a  better  estimate  of  the  surface  recession  and  surface  temperature. 
Both  analyses  do  not  account  for  possible  surface  gouging  caused  by  flow  sep¬ 
aration  as  mentioned  in  Section  5.3.1.  A  rough  estimate  of  the  total  gouging 
effect  was  given  in  Figure  5-10. 

3.3.3  PG  Washer  Backwall  Insulation  Effects 

To  estimate  the  effects  of  the  PG  backwall  boundary  condition  on  the 
thermal  response  of  the  PG  washers,  a  one-dimensional  analysis  was  performed 
using  an  insulated  PG  backwall  boundary  condition.  The  sensitivity  of  the  surface 
recession  to  this  boundary  condition  is  of  interest  due  to  the  possibility  of  a 
gap  between  the  washers  and  ATJ  backup.  Results  of  this  study  are  shown  in 
Figures  5-18  and  5-19.  Using  the  "new"  kinetics  and  the  thermochemistry  of  the 
18  percent  Al-CTPB  propellant  the  insulated  backwall  boundary  condition  caused  a 
small  increase  in  surface  temperature  and  approximately  a  25  percent  increase  in 
surface  recession. 
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Figure  5-1.  Throat  Package  Geometry  for  2-D  AiuOvsi 
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Figure  5-5.  Thermal  Properties  of  ATJ  Bulk  Graphite 
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Figure  5-6.  Estimated  Thermal  Properties  of  MX4926  Accounting 
For  Internal  Decomposition 
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Figure  5-9.  Continued 

B)  Exposure  Time  =  1.0  Seconds 
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Figure  5-9.  Continued 

C)  Exposure  Time  =  2.0  Seconds 
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Figure  5-12.  Comparison  of  Predicted  Surface  Temperatures  at  the 
Throat 
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igure  5-17.  Comparison  of  MX4926  Surface  Temperatures  Given  By 
The  one  and  Two  Dimensional  Analysis 
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SECTION  6 


PRETEST  NOZZLE  THERMAL  RESPONSE  PREDICTIONS 

One  of  the  objectives  of  contract  F04611-7.1-C-0051  is  to  demonstrate  the 
validity  of  existing  analysis  techniques  for  predicting  the  performance  of 
materials  used  in  high  chamber  pressure  nozzle  designs.  The  criteria  used  for 
satisfying  this  objective  was  to  make  pretest  predictions  at  18  locations  in 
selected  nozzles  of  the  initial  13  test  firings  which  were  tested  by  UTC  under 
under  the  above  mentioned  contract. 

Two  types  of  pretest  predictions  were  made.  From  those  locations  analyzed 
for  motors  S/N  003,  S/N  005,  and  S/N  006,  the  pretest  predictions  were  made 
using  the  predicted  chamber  pressure  history.  For  the  remainder  of  the  locations, 
the  pretest  predictions  were  made  using  the  post  test  chamber  pressure  histories. 
However,  the  reason  that  both  types  of  predictions  are  referred  to  as  "pretest"  is 
because  the  results  of  all  predictions  were  transmitted  to  UTC  prior  to  receiving 
the  post  test  measurements. 

The  surface  thermochemistry  maps  required  for  the  predictions  are  presented 
in  Section  6.1  for  all  the  materials  except  the  PG  washer.  A  summary  (except 
for  the  PG  washer)  of  the  predictions  and  their  comparison  to  post  test  measure¬ 
ments  of  surface  recession  and  char  depth  are  given  in  Section  6.2.  The  compari¬ 
son  of  the  prediction  and  measured  surface  recessions  for  the  PG  washer  is 
presented  in  Section  6.3  Finally  a  summary  of  the  post  test  char  density  measure¬ 
ments  is  presented  in  Section  6.4. 

6.1  SURFACE  THERMOCHEMISTRY  DATA 

The  thermochemical  calculations  performed  for  the  pretest  predictions 
are  identified  in  Table  6-1.  This  table  summarizes  the  locations,  materials, 
and  propellants  for  which  thermochemistry  maps  were  made  and  their  corresponding 
figures. 

The  majority  of  the  thermochemistry  calculations  are  for  charring  ablation 
materials.  For  transient  analyses,  these  materials  require  themochemistry  maps 
with  parametrically  varying  values  of  as  shown  in  the  figures.  A  transient 
map  was  not  made  for  R-155  since  sufficient  chemical  data  was  not  available. 

A  steady  state  map  is  given  in  Figure  4-1  for  R-155  and  was  used  to  generate  a 
series  of  steady  state  response  solutions  using  the  SSA  computer  code.  Collectively, 
these  solutions  provided  a  pseudo-transient  prediction  which  was  compared  to 
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to  the  measured  test  data.  Thermochemistry  maps  for  the  two  non-decomposing 
materials,  G-90  (5  percent  and  18  percent  AL-CTPB)  and  Pycobond,  are  also 
given  in  these  figures.  The  maps  for  G-90  were  generated  using  the  GASKET 
computer  code  which  incorporates  reaction  kinetics  in  the  surface  thermochemistry. 
(Reference  2).  Figure  6-9  shows  the  surface  response  of  Pycobond  to  5  percent 
Al-CTPB.  Reaction  kinetics  were  not  incorporated  in  this  map  since  previous 
studies  (Reference  1)  have  shown  that  the  surface  recession  of  carbon/carbon 
composites  can  be  predicted  with  an  accuracy  of  +  25  percent  using  diffusion 
controlled  thermochemistry. 


6.2  COMPARISON  OF  PREDICTED  AND  MEASURED  SURFACE  RECESSION,  CHAR  DEPTHS  AND 
MATERIAL  DENSITY  PROFILES 


This  section  summarizes  the  18  pretest  predictions  made  for  a  total  of 
six  test  firings.  All  but  one  prediction  was  made  using  the  CMA  computer  code. 
The  prediction  for  R-155  was  made  by  a  series  of  steady  state  predictions  using 
the  SSA  computer  code.  Input  for  these  codes,  as  before,  was  derived  from  the 
information  presented  in  Section  3.  The  heat  transfer  coefficients  were  obtained 
from  the  figures  in  Section  3.5.1  and  were  modified  by  Equation  (3-3)  which 
is  repeated  here. 


°eueCH  =  °-75 
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peueCH 


ARGEIBL 
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The  instantaneous  chamber  pressure  (Pc)  for  each  test  are  presented  in  Figures 
6-13  through  6-18.  The  actual  chamber  pressure  histories  shown  for  test  firings 
SN  003,  SN  005,  and  SN  006  were  not  used  in  making  the  pre-test  predictions, 
but  rather  the  predicted  chamber  pressure  histories  (supplied  by  UTC)  were  used. 
These  are  shown  along  with  the  actual  tracings  in  Figures  6-14  through  6-16. 

The  pressure  history  for  test  firing  SN  006  (Figure.  6-16)  shows  that  at  8 
seconds  after  start-up  the  pressure  decreased  rapidly.  This  was  due  to  the 
failure  and  subsequent  loss  of  the  throat  section.  Since  predictions  were  made 
using  the  higher  predicted  pressure  history  they  were  not  compared  to  the 
measured  results.  Figure  6-18  illustrates  the  pressure  history  of  test  firing 
13.  Tnis  nozzle  had  a  keyhole  slotted  grain  propel'.ant  which  caused  the  pressure 
to  reach  a  maximum  immediately  after  startup  and  then  continuously  decrease 
during  the  test. 

Table  6-2  summarizes  the  post  test  measurements  made  by  UTC  for  the 
locations  at  which  predictions  were  made.  Surface  recession  measurements, 
in  general,  were  made  in  more  than  one  circumferential  location  and  are 
averaged  to  provide  a  comparison  with  predictions.  In  the  blast  tube,  recession 
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measurements  were  also  averaged  over  two  different  axial  locations.  Additional 
measurements  were  provided  within  0.25“  of  the  blast  tube  entrance  and  exit,  but 
were  discarded  because  of  undesirable  end  effects.  Char  depth  measurements 
were  net  obtained  at  different  circumferential  locations,  but,  like  recession 
data,  were  given  at  different  axial  locations  in  the  blast  tube.  Averages 
of  measurements  taken  at  axial  distances  of  3.0  and  5.0  inches  were  used  for 
comparisons  with  the  analytical  predictions. 

A  summary  of  the  pretest  predictions  and  their  comparisons  to  measured 
aata  is  given  in  Tabic  6-3.  For  the  six  test  firings  listed,  a  total  of  18 
predictions  were  made  for  the  various  locations  and  materials  shown.  The 
thermochemical  recession  predictions  shown  are  those  obtained  from  the  CMA 
computer  code  (SSA  computer  code  for  prediction  16)  and  do  not  account  for  char 
swell.  Reference  1  was  used  to  estimate  the  char  swell  for  MXCE280  and  MX4926 
(same  as  FM5055).  Char  swell  information  for  the  other  decomposing  materials 
was  not  available.  The  predicted  net  recession  shown  in  Table  6-3  equals  the 
predicted  thermochemical  recession  less  the  amount  of  char  swell  listed.  The 
predicted  char  depth  was  also  that  obtained  from  the  CMA  computer  program, 
with  the  amount  of  char  swell  added.  Average  measured  surface  recessions  are 
those  taken  from  Table  6-2  except  for  predictions  5,  10  and  11.  The  exit 
cone  measurements  shown  in  Table  6-3  are  for  an  area  ratio  of  5.0  and  are  inter¬ 
polated  values  between  those  listed  in  Table  6-2  C  for  area  ratios  of  4.2  and 
6.7.  As  mentioned  before,  the  measurements  for  test  firing  S/N  006  are  not 
included  in  this  table  due  to  nozzle  failure.  Test  firing  S/N  013  used  a 
slotted  grain  propellant  and,  as  shown  in  Figure  3-4,  the  predicted  heat  transfer 
coefficient  distribution  is  influenced  by  the  two  dimensional  flow  effects. 

For  this  reason,  the  prediction  is  dependent  on  whether  the  analysis  location 
is  in  or  opposite  to  the  plane  of  the  slot.  Also  listed  in  Table  6-3  are  the 
average  measured  char  depths  of  which  those  for  tests  S/N  006  and  S/N  0013  were 
not  available. 

The  last  column  of  Table  6-3  summarizes  the  comparisons  between  measured 
and  predicted  surface  recessions.  Excellent  agreement  with  data  was  found, 
with  most  predictions  falling  withing  a  +  25  percent  error  band  about  the 
measured  values.  Prediction  5  fell  below  the  measured  value  but  it  should  be 
noted  that  the  measured  surface  recession  was  a  relatively  small  number  which 
had  an  uncertainty  band  larger  than  the  average  value  reported.  Predictions  8 
and  9  would  fall  within  5  percent  and  28  percent,  respectively,  <f  the  extremes 
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of  their  measured  uncertainty  bands  were  considered.  The  prediction  techniques 
used  for  the  slot  plane  in  test  13  is  shorn  to  be  reasonably  accurate  while 
further  studies  may  be  required  for  modelling  the  flow  field  in  the  plane 
opposite  the  slot. 

Typical  surface  recession  rate  and  surface  temperature  histories  are 
shown  in  Figure  6-19  for  an  MXCE280  material  in  the  blast  tube.  As  shown  in 
the  figure,  the  surface  recession  rate  decreases  during  the  motor  firing,  but 
the  surface  temperature  remains  approximately  constant. 

6.3  PG  WASHER  SURFACE  RECESSION 

The  surface  recession  rates  for  the  PG  washer  nozzle  throat  are  summa¬ 
rized  in  Figure  6-20  for  both  the  motors  tested  under  the  current  program 
(Contract  F04611-71-C-0051)  which  used  a  CTPB  propellant  with  5  and  18  percent 
aiiiminum  loading  fractions  and  those  tested  under  the  previous  program  (Contract 
F-4611-69-C-0065)  which  used  a  PBAN  propellant  with  16  percent  aluminum  loading 
fraction.  The  difference  in  the  measured  surface  recession  rate  for  the  two 
types  of  propellant  is  obvious.  The  predicted  surface  recession  rates  for  botn 
types  of  propellant  are  also  tabulated  on  this  figure.  These  rates  are  based 
on  the  "new"  kinetics.  Good  agreement  is  obtained  between  the  mear,  ired  and 
predicted  surface  recession  rates  for  the  PBAN  propellant  but  poor  agreement 
is  obtained  for  the  CTPB  propellant.  Before  going  into  further  discussion  on 
this  discrepancy  a  comparison  of  the  thermochemical  properties  of  the  two 
types  of  propellant  might  prove  beneficial.  This  comparison  is  made  for  the 
16  percent  aluminum  PBAN  propellant  and  the  18  percent  aluai.  um  CTPB  propellant. 

The  elemental  constituents  of  the  two  propellants  are  summarized  in  Table 
6-4.  The  gas  only  elemental  constituents,  which  were  obtained  by  removing  the 
A1203  particles  from  the  combustion  products,  are  summarized  in  Table  6-5.  The 
particle  concentration  was  defined  by  performing  an  ACE  chemical  equilibrium 
solution  at  the  conditions  shown  m  Table  6-6.  The  propellant  ballistics 
quoted  by  the  manufacturer  are  summarized  in  Table  6-7.  The  remaining  thermo¬ 
chemical  information  for  the  two  propellants  is  summarized  in  Table  6-8.  This 
table  presents  the  main  gaseous  products  for  the  two  propellants.  These  were 
defined  for  the  conditions  summarized  in  Table  6-6  and  were  obtained  by  perform¬ 
ing  an  ACE  chemical  equilibrium  solution.  A  comparison  of  the  gaseous  products 
for  the  two  propellants  shows  a  maximum  difference  of  approximately  15  percent 
for  the  main  oxidizing  species.  As  a  result,  the  chemical  constituents  of  the 
two  propellants  do  not  provide  an  obvious  clue  to  the  discrepancy  between  the 
predicted  and  measured  surface  recession  rates.  Another  possibility  is  the 


difference  in  temperatures  of  the  two  propellants  but  this  is  taken  into 
account  in  the  thermochemical  predictions.  The  effect  on  the  surface  recession 
rate  resolutions  from  the  difference  in  the  flame  temperatures  is  approximately 
a  factor  of  two.  This  factor  is  based  on  a  change  in  temperature  of  approximately 
300°R  from  a  temperature  of  500D°R  {approximate  surface  temperature) .  However, 
the  predicted  surface  recession  rate  shown  in  Figure  6-20  is  less  for  the  18° 

Al-CTPB  propellant  than  the  16°  Al-PBAN  propellant.  This  is  due  to  the  predictions 
for  the  CTPB  propellant  being  used  for  a  10  second  motor  firing  compared  to  a 
20  second  motor  firing  for  the  predictions  for  the  PBAN  propellant.  The  longer 
motor  firing  time  results  in  a  higher  average  surface  temperature  which  in  turn 
yields  a  higher  average  B'  and  surface  recession. 

Another  possible  explanation  for  the  difference  is  that  a  step  occurs 
at  the  junction  between  the  leading  PG  washer  and  the  G-90  aft  entrance  cap. 

This  step  could  have  the  effect  of  restarting  the  boundary  layer.  A  surface 
recession  prediction  was  made  for  this  condition,  and  it  resulted  in  a  factor  of  250% 
increase  in  the  predicted  surface  recession  rate.  This  is  also  considerably 
less  than  the  measured  surface  recession  rate. 

From  'he  above  comments,  it  is  apparent  that  the  performance  of  the  PG 
washer  in  the  CTPB  propellant  environments  is  not  understood.  The  large  slope 
of  the  curves  shown  in  Figure  6-20  suggest  that  the  discrepancy  is  due  to  some 
phenomena  other  than  kinetically  controlled  surface  chemical  reactions. 

6-4  CHAR  DENSITY  PROFILE  MEASUREMENTS 

Additional  experi.- e  ital  data  which  can  be  used  to  evaluate  the  analytical 
techniques  is  the  post  test  char  density  profile.  Several  attempts  were  made  at 
making  these  measurements  for  the  MXCE  280  material  in  the  blast  tube  of  motoi 
S/N  009.  As  shown  in  Tabel  6-3,  the  predicted  char  depth  in  this  blast  tube  was 
only  .005  i"ches  wnich  makes  che  experimental  definition  of  the  char  density 
profile  extremely  difficult. 

The  shave  and  weigh  method  was  used  in  trying  to  make  the  char  density 
measurements.  This  method  involves  sectioning  a  sample  from  the  blast  tube, 
weighing  the  sample,  shaving  a  small  amount  from  the  surface,  and  then  reweighing 
the  sample.  This  process  was  repeated  until  the  sample  was  consumed.  The  sample 
used  for  the  measurements  was  taken  from  a  location  which  was  approximately  one 
inch  from  the  upstream  end  of  the  blast  tube.  The  sample  had  a  surface  area  of 
.234  x  487  inches  with  the  longer  length  being  parallel  to  the  blast  tube  center¬ 
line.  Because  of  the  small  char  depth,  only  a  small  amount  could  be  shaved  off  the 
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sample  during  a  single  step.  For  the  first  attempt,  a  series  of  5  mil  thick 
slices  was  removed  from  the  surface  of  the  sample.  The  reduced  char  densities 
showed  considerable  scatter  due  to  the  weight  removed  during  a  single  skaving 
operating  being  the  difference  of  two  relatively  large  numbers.  A  second 
attempt  was  made  with  the  thickness  of  the  slices  being  increased  to  10  mils. 
However,  the  reduced  char  densities  still  showed  considerable  scatter.  Further 
attempts  were  not  made  because  increasing  the  thickness  of  the  slices  would 
result  in  this  thickness  being  a  significant  fraction  of  the  char  depth. 


TABLE  6-1 

SUMMARY  OF  PRETEST  THERMOCHEMISTRY  MAPS 


Location 

Material 

Propellant 
(%  AL) 

Figure  No. 

Aft 

Durez  16771-1 

5 

6-1 

Closure 

Durez  16771-1 

18 

6-2 

MX2625 

5 

6-3 

R-155 

18 

4-1 

Aft 

Entrance- 

G-90 

5 

6-4 

Cap 

Blast 

MXCE  280 

5 

6-5 

Tube 

MACE  280 

18 

6-6 

MXSE  280 

5 

6-7 

G-90 

18 

4-3 

MX4926 

5 

6-8 

Pycobond 

5 

6-9 

Exit 

FM5055 

5 

6-10 

Cone 

FM5055 

18 

6-11 

MX2600 

5 

6-12 
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TABLE  6-2 


MEASURED  MATERIAL  RECESSION  AND  CHAR  DEPTHS 
A)  AFT  CLOSURE 


Circumfer- 

ential 

Measured 

Measurement 

Recession 

Location 

(in) 

Average 

Measured 

Recession 

(in) 


Measured 

Char 

Depth  (in) 


16771-1 


TOC 

TDC+45 


.142+. 016 


DUREZ 

16771-1 


TDC 

TDC+45 


.186+. 030 


MX2625 


TDC 

TDC+45 


.041 +.032 


TDC 

TDC+45 

TDC+180 

TOC+225 


.185 

.070 

-  .031 

-  .058 


a  average  not  considered  due  to  asymmetric  slotted  grain  effects 
b  measurements  not  available 


TABLE  6-2  (Continued) 


MEASURED  MATERIAL  RECESSION  AND  CHAR  DEPTHS 
B)  BLAST  TUBE 


W- 


I 


’ - 1 

Average 

Measured 

Average 

j 

Pleasured 

Measured 

Char 

Measured 

Test  ■  Area 

Measurement  Location 

Recession 

Recession 

Depth 

Char  Depth 

!  Firing  .  Ratio 

Mat‘1 

Axial(in) 

Circum(Deg) 

(in) 

(in) 

(in) 

(in) 

1  l 

2  2.0 

MXCE280 

3.0 

Tf\r 
l  iA» 

.384 

.145 

.140 

5.0 

TDC+45 

TDC 

.347 

.3?6 

355*’- 029 
-.041 

.143+. 003 

TDC+45 

.314 

3  1.3 

MXCE280 

3.0 

TDC 

.440 

TDC+45 

.455 

.,.,+.014 

•456-.016 

.103 

TDC+180 

.459 

TDC+225 

.470 

G-90 

5.0 

TOC 

.193 

TDC+45 

.199 

196- 004 

’  -.003 

N/A 

TDC+180 

.200 

TDC+225 

.193 

5  1.3 

MXCE280 

3.0 

TDC+180 

.471 

.020 

.030 

5.0 

TDC+225 

TDC+180 

.565 

.526 

fim+-247 

.603-132 

.025+. 005 
- 

TDC+225 

.850 

MXCE280 

3.0 

TOC 

TDC+45 

.422 

.700 

.561 +.139 

.025 

MX4926 

5.0 

TOC 

TDC+45 

.588 

.733 

.661+. 073 

.085 

6  2.0 

PYCOBOND 

3.0 

TOC 

.200 

TDC+45 

.190 

T0C+180 

.209 

TDC+225 

.218 

ono+.Q] 5 

N/A 

5.0 

TOC 

.197 

•203-.013 

TDC+45 

.199 

TDC+180 

.201 

T0C+225 

.211 

9  2.0 

MXCE280 

3.0 

TOC 

.356 

TDC+45 

.376 

.110 

TDC+180 

.347 

5.0 

TDC+225 

TDC 

.324 

.346 

.350+. 026 

.105+. 005 

TDC+45 

.375 

.100 

TDC+180 

.347 

TOC+225 

.325 

13  1.3 

MXCE280 

3.0 

TDC 

.491 

TDC+45 

.525 

TDC+180 

-.035 

TDC+225 

.441 

a 

5.0 

TDC 

.500 

T0C+45 

.514 

i 

TDC+180 

-.021 

TDC+225 

.443 

_ 

L 

a  average  not  considered  due  to  asymmetric  slotted  grain  effects 
b  distance  from  entrance  to  bUst  tube 
c  measurements  not  available 
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TABLE  6-2  (Concluded) 
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MEASURED  MATERIAL  RECESSION  AND  CHAR  DEPTHS 
C)  EXIT  CONE 


Test 

firing 

Mat’  1 

Area 

Ratio 

Circumfer¬ 

ential 

Measurement 

Location 

Measured 

Recession 

(in) 

Average 

Measured 

Recession 

(in) 

Measured 

Char 

Depth  (in) 

3 

FM5055 

4.2 

TOC 

.223 

TOC+45 

.216 

.148+. 075 

.120 

TDC+180 

.070 

CO 

t-x. 

o 

1 

TDC+225 

.082 

6.8 

TOC 

-.003 

TDC+45 

r-*. 

o 

o 

1 

-.005+. 002 

.150 

TDC+180 

1 

o 

o 

cn 

TOC+225 

-.005 

5 

FM5055 

4.2 

TDC 

.122 

.123+. 001 

.065 

TOC+45 

.124 

6.7 

TDC 

.025 

.039+. 014 

.080 

TDC+45 

.053 

MX2600 

4.2 

T0C+180 

.114 

.108+. 006 

.060 

TOC+225 

.102 

6.7 

TDC+180 

.032 

.032+. 000 

.080 

TDC+225 

.032 

S 


TABLE  6-4 

ELEMENTAL  PROPELLANT  CONSTITUENTS 


Element 

Concentration 
(atomic  wt/unit  mass) 

CTPB 

18%  Al 

PBAN 

16%  Al 

H 

.03684 

.03747 

C 

.00824 

.01030 

N 

.00602 

.00596 

0 

.02422 

.02436 

A1 

.00663 

.00593 

Cl 

_ 

.00601 

.00579 

TABLE  6-5 


ELEMENTAL  PROPELLANT  CONSTITUENTS 
GAS  ONLY 


Concentration 
(atomic  wt/unit  mass) 


CTPB 

18%  A1 

PBAN 

16%  A1 

.05359 

.05217 

.01199 

.01434 

.00876 

.00831 

.02186 

.02237 

.0072 

.00056 

.00875 

.00806 

TABLE  6-6 

CONDITIONS  FOR 
PROPELLANT  EVALUATIONS 


Propellant 

Pressure 

Temperature 

(Atm) 

(°K) 

18%  Ai-CTPB 

190.5 

3692 

16%  Al-CTPB 

170.0 

3485 

TABLE  6-7 

PROPELLANT  BALLISTICS 
QUOTED  BY  MANUFACTURER 


Propellant 

Pressure 

Flame 

Oxidation 

(psia) 

Temperature 

(°K) 

Ratio 

18%  Al-CTPB 

1000 

3594 

1.33 

16%  Al-CTPB 

1000 

3420 

1.27 

TABLE  6-8 
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PROPELLANT  CONSTITUENTS 
(MOLE  FRACTION) 


Species  a 

181  Al-CTPB 

10%  Al-PBAN 

A1 

•160  x  10_3 

A1C1 

.526  x  10"3 

•364  x  1C"3 

aici2 

.130  x  10  1 

•111  x  10~* 

aiho2 

.110  x  10~3 

.487  x  10-3 

aici3 

.368  x  10"3 

.362  x  10" 3 

AlO 

.121  x  10"3 

.415  x  lO"1* 

CHO 

.346  x  10-3 

.340  x  10~3 

CO 

.226  x  10° 

.272  x  10° 

c°2 

.173  x  10“‘ 

.164  x  10"1 

Cl 

.110  x  10"' 

.692  x  10~3 

ClH 

.134  x  10° 

.129  x  10° 

HO 

.115  x  10~‘ 

.552  x  10~3 

H 

.340  x  10"' 

.239  x  10"3 

HP 

.279  x  10° 

.301  x  10° 

h2° 

.176  x  10° 

.145  x  10° 

NO 

,994  x  10_3 

.402  x  1 0~3 

N2 

.885  x  10-’ 

.834  x  10"1 

0 

.785  x  10~3 

.250  x  JO"3 

°2 

.195  1*T3 

— 

Al2°3 

.801  x  10_1 

.754  x  10" 1 

Includes  only  those  species  with  mole  fractions 
equal  to  or  greater  than  ,1  x  10“3 
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D I MENS 


Figure  6-2.  Surface  Response  of  Durez  16771  and  CTPB  18% 
Propellant 


NON-DIMENSION 


Figure  6-4.  Kinetically  Controlled  Response  of  G-90  Bulk  Graphite  to 
5%  AL-CTPB  Propellant 


rigure  6-5.  Surface  Response  of  MXCE  280  Carbon  Phenolic  to 
5«  AL-CTPB  Propellant 
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Figure  6-11.  Surface  Response  of  FM5055  to  CTPB  18% 
Propellant 


TEMPERATURE,  *R 

Figure  6-12.  Surface  Response  of  MX2600  Silica  Phenolic 
to  ?i%  AL-CTPB  Propellant 
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SECTION  7 


EXPERIMENTAL  EVALUATION  OF  THE  HEAT  FLUX 
BOUNDARY  CONDITIONS 

Total  heat  flux  sensors  (or  thermocoupled  plugs)  were  installed  in  the 
blast  tube  of  rocket  motors  S/N  005  and  S/N  007.  The  objectives  of  this  instru¬ 
mentation  were  to 

•  Determine  experimentally  the  total  heat  flux  to  the  nozzle  wall  and 
relate  this  to  the  convective  heat  transfer  coefficient  predicted 
by  the  anlaytical  model. 

•  Determine  experimentally  the  transient  surface  recession. 

The  thermal  instrumentation  systems  used  to  satisfy  the  above  objectivves, 
and  then  installation  into  the  nozzle  components  are  described  briefly  in 
Section  7.1.  The  data  reduction  procedure  is  summarized  in  Section  7.2,  and 
the  experimental  data  and  comparisons  between  the  data  and  analytical  predictions 
are  presented  in  Section  7.3. 

7.1  THERMAL  INSTRUMENTATION  SYSTEMS 

The  thermal  instrumentation  on  the  liyh  chamber  pressure  nozzle  firings 
consisted  of  four  total  heat  flux  sensors  and  three  narrow  view  angle  radio¬ 
meters.  The  instrumentation  was  concentrated  on  motors  S/N  0 05  and  S/N  007. 

This  instrumentation  is  described  in  considerable  detail  in  References  1,  19, 
and  20  and,  thus,  .ts  main  features  will  only  be  summarized  in  the  following 
paragraphs.  A  description  of  the  total  heat  flux  sensor  and  its  installation 
into  the  nozzle  components  will  be  presented  in  the  following  paragraphs. 

A  schematic  of  the  total  heat  flux  sensor  is  shown  in  Figure  7-1.  This 
sensor  is  0.500  inches  in  diameter  and  is  fabricated  of  the  same  material  at 
the  same  layup  angle  as  in  the  nozzle  component  in  which  it  is  to  be  installed. 
Thermocouple  lead  wires  are  placed  m  an  isothermal  plane  at  various  depths 
below  the  surface  with  the  thermocouple  junction  at  the  center  of  the  sensor. 

This  junction  which  is  nominally  15  mils  in  diameter  butts  up  against  the  sen¬ 
sor  material  at  the  interface  between  the  10  and  20  mil  cylindrical  holes.  Tne 
thermocouple  leads  are  placed  along  the  sides  of  the  sensor  and  are  electrically 
insulated  from  the  sensor  by  beryl lia  sleevinq  which  is  potted  in  a  boron  nitride 
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slurry.  The  leads  are  attached  to  30  mil  diameter  pins  located  at  the  back  of 
the  sensor.  These  pins  provide  terminals  for  attaching  the  measurement  system 
leads.  The  thermocouples  are  either  tungsten  5  percent  rhenium/tungsten  26 
percent  rhenium  or  Chrome 1 /A lumel  and  are  accurately  located  using  x-ray  photo- 
araphs. 

The  sensors  are  installed  in  the  ablative  nozzle  components  prior  to 
the  assembly  of  these  components  into  the  structural  steel  shell.  The  sensors 
are  bonded  to  the  nozzle  component  along  the  cylindrical  surface  of  the  sensor, 
and  in  addition,  the  cavity  behind  the  sensor  is  filled  with  epoxy  to  protect 
the  pins  and  the  thermocouple  lead  wire  junctions. 

The  locations  of  the  total  heat  flux  sensors  in  the  blast  tube  of  the 
two  nozzles  are  shown  in  Figure  7-2.  The  sensors  are  at  the  same  axial  location 
in  the  two  blast  tubes  but  their  radial  locations  were  different  due  to  the 
difference  in  the  initial  blast  tube  area  ratio  for  the  two  nozzles. 

7.2  DATA  REDUCTION  PROCEDURE 

The  data  reduction  procedure  for  defining  the  total  neat  flux  to  the  ex¬ 
posed  surface  of  a  solid  rocket  nozzle  is  outlined  m  Figure  2-3.  In  this 
procedure,  the  ACE  computer  program  is  used  to  define  the  wall  enthalpies  and 
the  non-dimensional  surface  thermochemical  maps,  and  the  Option  2  version  of  the 
CMA  computer  program  is  used  to  determine  the  in-depth  temperatures  and  the 
pertinent  in-depth  energy  terms.  The  in-depth  energy  terms  of  interest  are 
the  thermal  conduction,  the  adsorption  of  energy  by  the  pyrolysis  gases  as 
they  percholate  to  the  surface,  and  the  adsorption  of  energy  in  the  pyrolysis 
zone  due  to  material  degradation.  In  the  Option  2  version  of  the  CMA  program, 
the  surface  boundary  conditions  consist  of  a  specified  surface  temperature  and 
a  specified  surface  recession  rate.  As  shown  in  Figure  7-3,  these  surface 
boundary  conditions  are  varied  until  agreement  is  obtained  between  the  measured 
and  predicted  in-depth  temperatures.  In  addition  the  integral  of  the  specified 
surface  recession  rate  must  be  in  agreement  with  the  measured  surface  recession. 
The  surface  and  in-depth  temperatures  and  the  surface  recession  rate  defined 
from  the  final  iteration  using  the  Option  ?  version  of  the  CMA  computer  code 
are  referred  to  as  "determined"  parameters. 

As  shown  in  Figure  3.1,  the  input  to  the  CMA  and  ACE  computer  programs 
require  physical,  thermal,  and  chemical  material  property  data.  Since  this 
information  is  used  in  determining  the  total  heat  flux  to  the  surface,  the  use 
of  the  thermocoupl ed  ablative  plug  as  a  total  heat  flux  sensor  requires  that 
the  ablative  material  be  characterized.  The  material  used  to  fabricate  the 
sensors  for  the  rocket  motor  firings  was  MXCE  280.  The  thermal,  physical, 
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and  chemical  properties  for  this  material  is  presented  in  Reference  1.  The  sur¬ 
face  thermochemical  maps  and  wall  and  free  stream  enthalpies  which  were  generated 
using  the  ACE  computer  code  are  presented  in  Appendix  A.  The  use  of  these  curves 
in  determining  the  total  surface  heat  flux  is  described  in  detail  in  Reference 
19  and  will  not  be  repeated  here. 

7.3  COMPARISONS  BETWEEN  MEASURED  AND  PREDICTED  SURFACE  HEAT  FLUX 

Experimental  thermocouple  data  were  obtained  only  from  the  heat  flux  sen¬ 
sors  installed  in  the  blast  tube  for  S/N  005.  The  determined  in-depth  tempera¬ 
tures  are  compared  with  the  measured  temperatures  in  Figure  7-4  for  the  sensor 
in  the  forward  end  of  the  blast  tube  and  in  Ficure  7-5  for  the  sensor  installed 
in  the  af*-  end  of  the  blast  tube.  The  net  surface  recession  r<.  es  used  in  the 
data  ’.eduction  analysis  to  obtain  the  determined  in-depth  temperatures  are 
presented  in  Figure  7-6  for  the  forward  plug  (or  sensor) .  T-ve  shape  of  this 
curve  is  due  to  surface  recession  in  the  blast  tube  causing  the  incident  surface 
heat  flux  to  decrease  during  the  motor  firing.  A  comparison  of  this  determined 
surface  recession  (option  2  CMA)  and  the  predicted  surface  recession  (option  1 
CMA)  is  shown  in  Figure  7-7. 

The  end  result  of  the  data  rod’.ction  is  the  determined  convective  heat 
transfer  coefficient  which  is  compared  with  the  predicted  value  in  Figure  7-8. 

The  comparisons  of  bath  the  surface  recession  and  the  heat  transfer  coefficient 
is  within  20  percent,  and  the  determined  values  are  generally  low.  However,  the 
good  agreement  which  is  obtained  substantiates  the  analytical  techniques  being 
used . 
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Figure  7-3. 


Procedure  for  Defining  Total  Heat  Fluxes 
Using  Measured  In-Depth  Temperature  Histories 
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SECTION  8 


CONCLUSIONS  AND  RECOMMENDATIONS 

•  Aerotherm  analysis  techniques  (ACE,  ARGEIBL,  OMA,  ASTHMA,  and  GASKET) 
have  been  validated  for  predicting  the  thermochemical  performance  of 
charring  ablating  and  bulk  graphite  materials. 

•  The  predicted  surface  recession  for  the  pyrolytic  graphite  washers  at 
the  nozzle  throat  was  low  by  a  factor  of  at  least  five  for  the  18 
percent  Al-CTPB  and  5  percent  Al-CTFB  propellants.  However,  the 
analytical  techniques  had  predicted  accurately  the  surface  recession 
of  this  material  for  a  16  percent  Al-PBAN  propellant  at  a  similar 
chamber  pressure.  The  exact  reason  for  the  discrepancy  has  not  been 
identified.  Because  of  the  severity  of  the  problem,  further  analytical 
and  experimental  efforts  should  be  pursued. 

•  The  predicted  circumferential  variation  in  the  blast  tube  surface 
recession  rate  was  not  in  good  agreement  with  the  experimental 
measurements  for  the  keyhole  port  grain  configuration.  However,  this 
discrepancy  is  due  to  the  assumption  made  in  the  flow  :  .eld  analysis 
that  the  two-dimensional  flow  effects  resulting  from  the  keyhole  port 
would  be  damped  out  in  the  upstream  end  of  the  blast  tube.  The  ex¬ 
perimental  data  show  that  the  two  dimensional  flow  effects  exist  along 
the  center  length  of  the  blast  tube. 

•  The  total  heat  flux  sensors/or  thermocouple  plugs)  provide  a  reliable 
and  accurate  system  for  defining  surface  recession  rates  and  in-depth 
temperatures  as  well  as  total  incident  surface  heat  flux.  The  convective 
heat  transfer  coefficient  determined  from  these  measured  data  were  in 
good  agreement  with  the  predicted  values. 

•  Thermochemical  screening  provides  an  economical  and  practical  method  for 
selecting  candidate  materials. 
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EXPOSURE  TIME,  SECONDS 


APPENDIX  A 


THERMOCHEMICAL  MAPS  FOR 
DEFINING  TOTAL  oCJRFACE  HJAT  FLUX 

The  thermcchemical  maps  required  for  defining  the  total  surface  heat 
flux  from  measured  in-depth  temperatures  are  presented  in  Figures  A-l  through 
A-4  for  MXCE  280.  The  information  presented  in  these  figures  include  the 
following: 

•  Relationship  between  non-dimensional  char  removal  rate,  (B^_)  non- 
dimensional  pyrolysis  off  gas  rate,  (B^)  ratio  of  B^/B and  wall 
temperature 

•  Relationship  between  the  wall  chemical  potential,  the  wall  temper¬ 
ature  and  the  non-dimensional  pyrolysis  off  gas  rate. 

•  Relationship  between  the  wall  enthalpy,  the  wall  temperature,  and 
the  non-dimensional  pyrolysis  off  gas  rate. 

•  Relationship  between  the  freestream  chemical  potential  and  the 
freestream  enthalpy  with  temperature. 

•  Recovery  enthalpy. 
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Figure  a~2  *  Non-Dimensional  Pyrolysis  Off  Gas  Rate  as  a  Function  of  Mass  Flux 
Ratio  for  MSCE  280  Blast  Tube  Sensor 


Freestream  and  Wall  Chemical  Potential  as  a 
Function  of  Wall  Temperature  for  MXCE  280 
Blast  Tube  Sensor. 
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Figure  A_<j 
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WALL  TEMPERATURE,  °R 

Recovery,  Freestream,  and  Wall  Enthalpies  as  a 
Function  of  Wall  Temperature  for  MSCE  280  Blast 
Tube  Sensor 
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